Systems and methods for bi-static or multi-static holographic navigation

ABSTRACT

The application relates to bi-static or multi-static holographic navigation systems, including methods of localizing an emitter or receiver with high precision relative to the sea floor. The system and methods can be used with a fully active sonar or radar system using well synchronized transmitters and receivers. The system and methods can be used with a passive sonar or radar system localizing a transmitter or a receiver based on poorly timed received signals.

CROSS-REFERENCE TO RELATED APPLICATIONS

This application is a continuation of U.S. patent application Ser. No. 13/848,570, filed Mar. 21, 2013, and entitled “Systems and Methods for Bi-Static or Multi-Static Holographic Navigation,” which claims priority to and the benefit of U.S. Provisional Patent Application Ser. No. 61/613,838, filed on Mar. 21, 2012 and entitled “Bi-Static or Multi-Static Holographic Navigation,” the entire contents of each of which are incorporated herein by reference.

FIELD OF THE INVENTION

The present disclosure relates generally to systems and methods associated with synthetic aperture sonar (SAS) technology. More particularly, in various implementations, the present disclosure relates to holographic navigation.

BACKGROUND

In most land-based applications, navigation is often aided by in-place infrastructure such as GPS, radio beacons or a priori maps. Navigation and mapping underwater is difficult because among other things, wide-coverage underwater GPS-equivalents do not exist and large portions of the sea bed are still unexplored.

Current techniques for underwater navigation use publicly available bathymetry maps. However, these maps are relatively coarse and unsuitable for precision navigation. Other sonar-based navigation systems rely on positioning schemes that use the sonar data itself. For example, on-the-fly acoustic feature-based systems attempt to use sonar to detect naturally occurring landmarks. Other solutions to the navigation problem include deploying low-cost transponders in unknown locations thereby enabling range-based measurements between the vehicle and transponder beacon. However, these transponders are often deployed at locations that are at great distances from each other, and often only partially observable because of the range-only information. Thus, these technologies are unsuitable for navigation across small vehicle paths.

It is often desirable to be able to navigate terrain (whether on land or underwater) in a vehicle equipped appropriately with sensors that allow the vehicle to navigate the terrain relative to a prior map of the terrain. Holographic navigation is a technique that allows for precise navigation of a vehicle equipped with sonar sensor array(s) relative to a prior map. Quite often, these vehicles are autonomous underwater vehicles (AUVs) or unmanned aerial vehicles (UAVs).

A synthetic aperture sonar (or radar) array on such a vehicle generates at least one image that can be compared to an image associated with the map of the terrain and may be processed by a computer system on the vehicle to navigate the vehicle in relation to the prior map of the terrain. For instance, SAS arrays enable coherent correlation between sonar signals, whether generated by the synthetic aperture, or not.

Existing holographic navigation systems suffer from a plurality of deficiencies, including the amount of power consumed by, the size, and the shape of these systems. In some existing systems, for example, the performance of a holographic navigation algorithm degrades substantially as transmitter frequencies increase and wavelengths decrease.

Furthermore, accurately localizing a transmitter based on received emissions is a known problem. Direct path observations can be used with an array or antenna to get a bearing to an object. Multiple direct path observations can be used for triangulation while waveguide models can be used for ranging from a single location. However, these techniques are all limited in precision. Extremely long arrays (such as SOSUS) can focus on the location of an emitter object with fairly high precision provided that it is in the near field, but are extremely expensive. Hence, there is not an existing method of passively localizing an emitter with extremely high precision using a low cost receiver.

Accordingly, there is a need for high precision and low cost navigation systems, particularly for underwater applications.

SUMMARY

The devices, systems and methods of the present application address these and other deficiencies of existing holographic navigation systems. Holographic navigation is an inexpensive way of coherently localizing an active, monostatic (or quasi-monostatic) sonar or radar relative to a prior map of terrain. Instead of using a physical aperture on the vehicle, the terrain acts as the aperture, yielding resolution that is a function of the spotlight on the seafloor. Given that the spotlight size typically grows with range, the scattering from that spotlight can be used to localize a robot with range invariant precision.

The present application includes a system and method of localizing an emitter or receiver with high precision relative to the sea floor. It can be used as either a fully active system; using well synchronized transmitters and receivers or as a passive system; localizing a transmitter or a receiver based on poorly timed received signals.

Placing a low frequency projector onto a separate platform (other than a receive platform), such as, without limitation, a large diameter unmanned underwater vehicle, may allow, for example, a mine neutralizer to reduce its sonar to small number of receive elements, which may have a relatively low weight. This enables low frequency holographic navigation on a very small platform and, ultimately, a more robust and flexible mine neutralization solution. Similarly, submarines prefer not to use active sonar as it may reveal their location, but they have very high quality receivers. Hence, bi-static holographic navigation can allow submarines to localize themselves based on seafloor scattering of a sonar signal originating from a known source (either fixed or moving).

Bi-static holographic navigation can also be used to localize a transmitter using a known receiver. Submarines and stealth aircraft often possess sonars and radars with electronically scanned phased arrays that steer energy into narrowly defined areas to reduce the probability of intercept. Although their outbound beams are narrow, the scattering off the seafloor or ground is less directional. Bi-static holographic navigation can detect and localize those target based on scattering. Similarly, ships and whales are difficult to accurately track and localize. However, bi-static holographic navigation can be used to track them and reduce collisions.

Bi-static holographic navigation, in its simplest form, may include two platforms with accurate synchronized clocks, a known transmit signal with known timing, position information about the transmitter, and a prior bi-static map. The receive system (either a single channel or an array) receives the signal, forms an image using timing and geometry information, applies grazing angle compensation, correlates the image against a prior image, converts the correlation result to a probability density function, and uses that probability density function to update its estimate of the state (position). High Frequency (HF) holographic navigation techniques can be applied to remove various forms of phase error.

In a more advanced form, the emitting system (submarine, aircraft, and the like) can emit a known signal from an unknown position at an unknown time that would be received at a known location at a known time. That signal would be back-projected or beam-formed, correlated, and used to update a position estimate.

In a more difficult instance, the signal being emitted by the submarine or aircraft would be unknown. In this case, if the direct arrival were known, it could be used to form an image. If the signal were unknown, it would be necessary to learn the signal, use it to form an image, and localize the target.

These techniques can also be performed multi-statically; applying equally to sonar and radar. Although back-scatter is more desirable, forward scatter can be used. Bi-static holographic navigation follows a similar process to monostatic holographic navigation. First, it is necessary to create a prior map. Ideally, this is a bi-static synthetic aperture image containing all possible frequencies at all possible angles. In practice, this will likely be band limited and angle limited. The remaining process depends on whether the transmitter and receiver are synchronized, which system has a known position, and whether the signal is known.

The system and methods herein have many applications. For instance, consider the case of localizing a vehicle based on transmissions from a known source. This applies to scenarios such as mine neutralization and submarine navigation. Bi-static holographic navigation allows small autonomous mine neutralizers to use low frequency signals. Low frequency acoustic signals are more robust to sedimentation than high frequency signals, but low frequency projectors can weigh more than neutralization platforms. Hence, by utilizing a separate transmitter platform having sufficient size to support a low frequency transmitter system, smaller neutralization platforms including receivers can be employed that navigate with high precision.

In one aspect, a system for navigating an underwater terrain includes a first platform. The first platform includes an acoustic transmitter arranged to transmit an acoustic signal toward a portion of the underwater terrain, resulting in a scattering signal emanating from the portion of the terrain. The first platform includes a first clock arranged to maintain timing information. The first platform also includes a memory arranged to store: i) position information associated with at least one of the acoustic transmitter and an acoustic receiver and ii) one or more settings associated with the acoustic signal. The first platform further includes a processor arranged to control the acoustic transmitter

The system also includes a second platform having an acoustic receiver arranged to receive a portion of the scattering signal. The second platform includes second clock that is synchronized with the first clock. The second platform also includes a memory arranged to store: i) position information associated with at least one of the acoustic transmitter and acoustic receiver, ii) the one or more settings associated with the acoustic signal, and iii) a monostatic or bi-static map of the underwater terrain where the map includes acoustic data, within a first frequency range, obtained from synthetic aperture sonar imaging of the underwater terrain. The second platform includes a second processor arranged to: i) generate a real aperture image of the portion of the underwater terrain based on the received portion of the scattering signal within a second frequency range and the timing information where the second frequency range can at least partially overlap with the first frequency range, ii) modify the real aperture image by compensating for grazing angle errors to generate a grazing angle invariant real aperture image and correct for phase errors in the grazing angle invariant real aperture image, iii) coherently correlate the modified real aperture image with the map, iv) convert the correlation result into a probability density function, and v) estimate the position of one of the first platform and the second platform based on the probability density function.

In one configuration, the first and second frequency ranges include a maximum frequency less than or equal to about 100 kHz. In other configuration, the first and second frequency ranges include a minimum frequency greater than about 100 kHz. The map may include a monostatic or bi-static map. Over modest angles, a monostatic prior map can be a substitute for a bistatic prior map. The first platform may be substantially stationary while the second platform is moving. The first platform may be moving while the second platform is substantially stationary. An acoustic signal may include one or more characteristics or settings such as, without limitation, a noise sequence, linear chirp, encoding, power, frequency, phase, and the like.

In some implementations, generating a real aperture image includes a plurality of real aperture images where each of the plurality of real aperture images represents a subset of the portion of the underwater terrain. The phase error in each of the plurality of real aperture images may be substantially constant. The acoustic receiver may include a single channel receiver or an array of receivers At least one of the first and second platforms may include an autonomous underwater vehicle (AUV). The underwater terrain may include at least a portion of a sea bed.

In another aspect, an underwater vehicle includes an acoustic receiver arranged to receive a portion of a scattering signal where the scattering signal emanates from a portion of an underwater terrain. The scattering signal may result from an acoustic signal being directed toward a portion of the underwater terrain by an acoustic transmitter on a transmitter vehicle. The vehicle also includes a first clock having timing information and being synchronized with a second clock located at the transmitter vehicle. The vehicle further includes a memory arranged to store: i) position information associated with at least one of the transmitting vehicle and the underwater vehicle, ii) one or more settings associated with the acoustic signal, and iii) a map (monostatic and/or bi-static) of the underwater terrain where the map includes acoustic data, within a first frequency range, obtained from synthetic aperture sonar imaging of the underwater terrain. The vehicle may also include a processor arranged to: i) generate a real aperture image of the portion of the underwater terrain based on the received portion of the scattering signal within a second frequency range and the timing information where the second frequency range at least partially overlaps with the first frequency range, ii) modify the real aperture image by compensating for grazing angle errors to generate a grazing angle invariant real aperture image and correcting for phase errors in the grazing angle invariant real aperture image, iii) coherently correlate the modified real aperture image with the map, iv) convert the correlation result into a probability density function, and v) estimate the position of one of the underwater vehicle and the transmitter vehicle based on the probability density function.

In some implementations, the first and second frequency ranges include a maximum frequency less than or equal to about 100 kHz. In other implementations, the first and second frequency ranges include a minimum frequency greater than about 100 kHz. The transmitter vehicle may be substantially stationary while the underwater vehicle is moving. The transmitter vehicle may be moving while the underwater vehicle is substantially stationary. Both vehicles may be moving or stationary.

In yet another aspect, a transmitter vehicle that enables underwater navigation by an underwater vehicle includes an acoustic transmitter arranged to transmit an acoustic signal toward a portion of an underwater terrain. This results in a scattering signal emanating from the portion of the terrain and a portion of the scattering signal being received by the underwater vehicle. The underwater vehicle uses the portion of the scattering signal to: i) generate a real aperture image of the portion of the underwater terrain based on the received portion of the scattering signal within a second frequency range and the timing information where the second frequency range at least partially overlapping with the first frequency range, ii) modify the real aperture image by compensating for grazing angle errors to generate a grazing angle invariant real aperture image and correct for phase errors in the grazing angle invariant real aperture image, iii) coherently correlate the modified real aperture image with the map, iv) convert the correlation result into a probability density function, and v) estimate the position of one of the underwater vehicle and the transmitter vehicle based on the probability density function.

The transmitter vehicle also includes a first clock arranged to maintain timing information that is synchronize with a second clock located at the underwater vehicle. The transmitter vehicle further includes a memory arranged to store: i) position information associated with at least one of the transmitting vehicle and the underwater vehicle and ii) one or more the settings associated with the acoustic signal. The transmitter also includes a processor arranged to control the acoustic transmitter.

The first and second frequency ranges may include a maximum frequency less than or equal to about 100 kHz. The first and second frequency ranges may include a minimum frequency greater than about 100 kHz. The transmitter vehicle may be substantially stationary while the underwater vehicle is moving. The transmitter vehicle may be moving while the underwater vehicle is substantially stationary. Both vehicles may be moving or stationary.

BRIEF DESCRIPTION OF THE DRAWINGS

The foregoing and other objects, features, advantages, and illustrative implementations of the invention will now be described with reference to drawings in which like reference designations refer to the same parts throughout the different views. These drawings are not necessarily to scale, emphasis instead being placed upon illustrating principles of the implementations.

FIG. 1 is a block diagram depicting a sonar mapping and navigation system.

FIG. 2 is block diagram of an exemplary computer system for implementing at least a portion of the systems and methods described in the present disclosure.

FIG. 3 depicts a transducer array in a sonar system.

FIG. 4 depicts a transducer array in a synthetic aperture sonar (SAS) system.

FIGS. 5A-5B depict a process for navigating a terrain using an exemplary high-frequency sonar navigation system.

FIGS. 6A-6B depicts a process for navigating a terrain using an exemplary high-frequency sonar navigation system.

FIG. 7 depicts a process for correcting range varying phase errors in a high-frequency sonar system.

FIG. 8 depicts a process for using a plurality of orthogonal signals in a synthetic aperture sonar (SAS) system to generate images.

FIGS. 9A and 9B depict a transducer array used in connection with an implementation of the process depicted in FIG. 8.

FIG. 10 depicts a process for transmitting pulses from a synthetic aperture sonar (SAS) system having multiple transmitters.

FIGS. 11A-C depict a transducer array used in connection with an implemenation of the process depicted in FIG. 10.

FIG. 12 depicts a process for simultaneous localization and mapping (SLAM) using real aperture sonar images.

FIG. 13 shows a bi-static holographic system including a transmitter and receivers operating in an underwater environment.

FIG. 14 is a flow diagram of a process for estimating a receiver position where timing is known, transmitter position is known, the acoustic signal is known, but the receiver position is unknown.

FIG. 15 is a flow diagram of a process for estimating a receiver position where timing is unknown, transmitter position is known, the acoustic signal is known and the receiver position is unknown.

FIG. 16 is a flow diagram of a process for estimating a receiver position where timing is unknown, transmitter position is known, the acoustic signal is unknown and the receiver position is unknown.

FIG. 17 is a flow diagram of a process for estimating a transmitter position where timing is known, transmitter position is unknown, the acoustic signal is known and the receiver position is known.

FIG. 18 is a flow diagram of a process for estimating a transmitter position where timing is unknown, transmitter position is unknown, the acoustic signal is known and the receiver position is known.

FIG. 19 is a flow diagram of a process for estimating a transmitter position where timing is unknown, transmitter position is unknown, the acoustic signal is unknown, but the receiver position is known.

DETAILED DESCRIPTION

To provide an overall understanding of the systems and methods described herein, certain illustrative implementations will now be described, including systems and methods for mapping and navigating a terrain. However, it will be understood by one of ordinary skill in the art that the systems and methods described herein may be adapted and modified for other suitable applications and that such other additions and modifications will not depart from the scope thereof.

The present application includes a system and method of localizing an emitter or receiver with high precision relative to the sea floor. It can be used as either a fully active system; using well synchronized transmitters and receivers or as a passive system; localizing a transmitter or a receiver based on poorly timed received signals.

Placing a low frequency projector onto a separate platform (other than a receive platform), such as, without limitation, a large diameter unmanned underwater vehicle, may allow, for example, a mine neutralizer to reduce its sonar to a small number of receive elements, which may have a relatively low weight. This enables low frequency holographic navigation on a very small platform and, ultimately, a more robust and flexible mine neutralization solution. Similarly, submarines prefer not to use active sonar as it may reveal their location, but they have very high quality receivers. Hence, bi-static holographic navigation can allow submarines to localize themselves based on seafloor scattering of a sonar signal originating from a known source (either fixed or moving).

The systems and methods described herein also include high-frequency (“HF”) holographic navigation, namely map-based navigation using the multi-aspect holographic-nature of synthetic aperture sonar (SAS) images captured at frequencies greater than or equal to about 100 kHz. The systems and methods described herein further include low frequency (“LF”) holographic navigation at frequencies less than about 100 kHz. In particular, the systems and methods described herein allow for coherent correlation between images, currently captured, and prior maps when there is an overlap in frequency and aspect. Such coherent correlation allows for position and/or heading-based navigation. At high-frequency, the inventor has recognized that images suffer from spatially varying phase errors (e.g., range varying phase errors), which cause image and/or correlation distortion. Such phase errors may exist even at low frequencies when there are altitude variations. In certain implementations, when the phase errors are much smaller than the bandwidth, although images may not be distorted, correlation (and therefore navigation) may become difficult. The systems and methods described herein overcome the deficiencies of the prior art by introducing a phase error corrector configured to cut the image into smaller regions where phase is relatively constant and use these phase measurements to correct portions of the image.

The systems and methods described here make use of various other aspects of the holographic nature of synthetic aperture images, which the inventor has recognized. For example, systems and methods are described herein for determining a three-dimensional model of a shape based on its two dimensional shading and shadowing of acoustic signals. The systems and methods described herein include methods for positioning sensors (such as Tsunami sensors) and navigation beacons with high-precision using HF holographic navigation. The systems and methods described herein include methods for monitoring and modeling a water column using an autonomous underwater vehicle (AUV) based on high-precision location measurements obtained using HF holographic navigation. In certain implementations, the systems and methods include a seismic survey system having a combination of orthogonal transmitters and multiple receivers to form a full planar synthetic aperture sonar with higher resolution.

In other aspects, the systems and methods described herein include adding multiple transmitters to the array and generating orthogonal pinging sequences. In particular, the systems and methods described herein include a SAS having a low-grating sidelobe, a SAS having a high coverage rate using multiple transmitters, and an overpinging sequence for increasing the range of the SAS system. The systems and methods described herein further include bistatic and monostatic holographic gapfilling techniques for localizing an emitter or receiver with high precision relative to a terrain. In still other aspects, the systems and methods described herein include simultaneous localization and mapping (SLAM) techniques that involve beamforming a real aperture image such that it can be coherently correlated with a prior real aperture image of overlapping frequencies. Each of these and other systems and methods described herein may be used independently of each other or in any suitable combination of one or more any other system and method. Modifications and variations described with reference to a system and method described herein may be applied to any other system and method described herein, without departing from the scope of the present disclosure.

In the following passages, an illustrative mapping and navigation system and an illustrative computer system for executing holographic navigation and mapping is described with reference to FIG. 1-4, respectively. Further illustrative implementations of components and processes of the holographic navigation and mapping system include processes for navigating a terrain, for example an underwater terrain, using a map are described with reference to FIGS. 5 and 6. To allow for high-frequency holographic navigation, FIG. 7 describes a process for correcting range varying phase errors, recognized by the inventor to be a reason for the failure of traditional holographic navigation, coherent correlation, and change detection systems at higher frequencies. FIGS. 8-9B depict a process and components for generating SAS images having low grating sidelobes, and FIG. 10-11B depict a process and components for generating a high-coverage rate SAS signals. Finally, FIG. 12 describes a holographic SLAM process for navigating a terrain.

FIG. 1 is a block diagram depicting a sonar mapping and navigation system 100. The system 100 includes a sonar unit 110 for sending and receiving sonar signals, a preprocessor 120 for conditioning a received (or reflected) signal, and a matched filter 130 for performing pulse compression and beamforming. The system 100 is configured to allow for navigating using high-frequency (greater than about 100 kHz) sonar signals. To allow for such HF navigation, the system 100 includes a signal corrector 140 for compensating for grazing angle error and for correcting phase error. The system 100 also includes a signal detector 150 for coherently correlating a received image with a map. In certain implementations, the system may be mounted on vehicle navigating over a terrain, such as an autonomous underwater vehicle (AUV) or an unmanned aerial vehicle (UAV). In such implementations, the system 100 includes an on-board navigation controller 170, motor controller 180 and sensor controller 190. The navigation controller 170 may be configured to receive navigational parameters from a GPS/RF link 172 (when available), an accelerometer 174, a gyroscope, and a compass 176. The motor controller 180 may be configured to control a plurality of motors 182, 184 and 186 for steering the vehicle. The sensor controller 190 may receive measurements from the battery monitor 172, a temperature sensor 194 and a pressure sensor 196. The system 100 further includes a central control unit (CCU) 160 that may serve as a hub for determining navigational parameters based on sonar measurements and other navigational and sensor parameters, and for controlling the movement of the vehicle.

In the context of a surface or underwater vehicle, the CCU 160 may determine navigational parameters such as position (latitude and longitude), velocity (in any direction), bearing, heading, acceleration and altitude. The CCU 160 may use these navigational parameters for controlling motion along the alongtrack direction (fore and aft), acrosstrack direction (port and starboard), and vertical direction (up and down). The CCU 160 may use these navigational parameters for controlling motion to yaw, pitch, roll or otherwise rotate the vehicle. During underwater operation, a vehicle such as an AUV may receive high-frequency real aperture sonar images or signals at sonar unit 110, which may then be processed, filtered, corrected, and correlated against a synthetic aperture sonar (SAS) map of the terrain. Using the correlation, the CCU may then determine the AUV's position, with high-precision and other navigational parameters to assist with navigating the terrain. The precision may be determined by the signal and spatial bandwidth of the SAS map and/or the acquired sonar image. In certain implementations, assuming there is at least a near perfect overlap of the sonar image with a prior SAS map with square pixels, and assuming that the reacquisition was performed with a single channel having a similar element size and bandwidth, and assuming little or no losses to grazing angle compensation, the envelope would be about one-half the element size. Consequently, in certain implementations, the peak of the envelope may be identified with high-precision, including down to the order of about 1/100^(th) of the wavelength. For example, the resolution may be less than 2.5 cm, or less than 1 cm or less than and about 0.1 mm in the range direction.

Generally, terrain recognition using long wavelength (low-frequency) sensors may be difficult due to the aspect dependence of object signatures. Sonar or radar images may be dominated by speckle that change with both sonar and object aspect, making incoherent image correlation extremely difficult. Coherently, any correlation operation involving signals with non-overlapping frequency bands will yield an answer of zero (since correlation is multiplication in the frequency domain). For two sonar images to correlate it is not enough that their spatial frequencies overlap, but the same points in the two images must be represented at overlapping frequencies. For a generic real aperture sonar, the same signature for a complex scene can only typically be re-observed by revisiting the original observation position and orientation and using the same frequencies. Consequently, in general, getting two complex sonar or radar images to coherently correlate is a measure zero occurrence; the expected cross correlation can be proven to be approaching zero. Therefore, coherently navigating relative to terrain is, in general, impossible if the system compares real aperture imagery to prior real aperture imagery, except as described below with reference to FIG. 12. Incoherent navigation is possible (i.e. using only the envelope) if there is distinct terrain, but against a uniform bottom (mud flat, field of gravel, ocean floor, etc.) this is usually not so.

Holographic navigation of a terrain, e.g., using a system implemented on AUVs, solves this problem by replacing at least one of the real aperture images with a synthetic aperture image. Because a synthetic aperture image is a type of hologram (or quasi-hologram) it contains all possible real aperture images over some range of frequencies and angles. Consequently, it may be possible to correlate a real aperture image against the synthetic aperture image and have a non-zero expected cross correlation. However, according to the Closed/Open Aperture theorem, it may be required that the synthetic aperture be a planar synthetic aperture, meaning that it is fully populated and Nyquist sampled in two dimensions. This type of population and sampling frequency is, in general, impractical.

By assuming the terrain is a manifold with embedded scatterers on the surface, and avoiding sub-bottom profiles/operating above the critical angle, or operating below the critical angle where the SNR is low, it is possible to show that the planar aperture can be replaced with a contour aperture provided the frequencies can rescaled. For example, consider an active sonar or radar and two scatterers spaced 5 centimeters apart in range on a flat bottom. From the perspective of a sonar or radar looking at the scatterers from the ground, the distance of travel for the two echoes differ by 10 cm (out and back). If the observer is, instead, looking down at an angle of 45 degrees above horizontal, the difference is shorted by cosine of 45 degrees (half) to 7.07 cm. So at horizontal a 10 cm wavelength would be exactly one cycle out of phase (constructively interferes), and a 20 centimeter wavelength would be exactly a half cycle out of phase (destructively interfere). At 45 degrees, the same would be true of a 7.07 cm wavelength and a 14.14 cm wavelength. Both wavelengths are scaled by the same amount (and, similarly, so are frequencies, except inversely). More generally, a change in vertical angle shifts all frequencies and changes the signal length by the cosine of the angle. This is not a shift in frequency so much as a change in pitch, where a doubling in frequency corresponds to a change in pitch of one octave. So by changing the observation angle from horizontal to looking down at 60 degrees the expected return is shorted by half and increases in pitch by one octave. In order for this to work, it is necessary for the second observation to be made with appropriately scaled frequencies relative to the first; for a very narrowband system too much of a change in grazing angle simply leads to the known signatures being out of band.

In some configurations, using grazing angle compensation and a prior synthetic aperture image of the systems and methods described herein, it is possible to navigate relative to terrain using a single element sonar or radar. Although synthetic aperture systems are extremely expensive, single element systems are generally very cheap. This means a very expensive mapping system can enable the widespread use of cheap autonomous systems with minimal inertial navigation. However, successful holographic navigation implementations to date have all used low frequency sonars (i.e. under 50 kHz), while the higher frequency systems have not worked. This is unfortunate, because lower frequency transmitters are, in general, larger, higher power, and more expensive. Thus, it is desirable to have a high frequency single element holographic navigation system. Further illustrative implementations of holographic navigation systems and methods are disclosed in U.S. patent application Ser. Nos. 12/802,453, 12/454,486, 12/454,484, 12/454,885, 13/466,059, 13/466,062, 13/466,063, 13/466,067, 13/466,067, 13/466/075, and 13/466,078, the contents of each of which are incorporated herein by reference in their entirety.

In one aspect, the invention relates to a method of terrain relative localization via holographic navigation. Holographic navigation and holographic maps are further described in U.S. patent application Ser. Nos. 12/798,169 and 12/802,455, the contents of each of which are incorporated herein by reference in their entirety. In some respects, holographic navigation is a method of terrain relative localization that takes advantage of the holographic properties of sonar and radar images. Quite often such terrain relative localization is performed by a system implemented on an autonomous underwater vehicle (AUV). However, the performance of holographic navigation algorithms implemented on such systems may degrade substantially as frequencies increase and wavelengths decrease. Conventionally, it is generally assumed that such degradation is because the some of the assumptions of grazing angle compensation break down. In other words, it is assumed that a change in vertical aspect no longer maps to a pure change in pitch because shadowing, occlusion, and complex three dimension relief fundamentally change the signature. However, the inventor has recently recognized that this assumption is not entirely incorrect, and that holographic navigation may fail at higher frequencies due to spatially varying phase errors. In some implementations, the invention corrects for those range varying phase errors by allowing for holographic navigation at higher frequencies with lower power consumption and smaller sized hardware.

As noted above, the system 100 includes a sonar unit 110 for transmitting and receiving acoustic signals. The sonar unit includes a transducer array 112 having a one or more transmitting elements or projectors and a plurality of receiving elements arranged in a row. In certain implementation, the transducer array 112 includes separate projectors and receivers. The transducer array 112 may be configured to operate in SAS mode (either stripmap or spotlight mode) or in a real aperture mode. In certain implementations, the transducer array 112 is configured to operate as a multibeam echo sounder, sidescan sonar or sectorscan sonar. One example of a transducer array is shown in FIG. 3 having one transmitting elements and six receiving elements. The transmitting elements and receiving elements may be sized and shaped as desired and may be arranged in any configuration, and with any spacing as desired without departing from the scope of the present disclosure. As described later in the present disclosure the number, size, arrangement and operation of the transducer array 112 may be selected and controlled to insonify terrain and generate high-resolution images of a terrain or object. One example of an array 112 includes a 16 channel array with 5 cm elements mounted in a 12¾ inch vehicle.

The sonar unit 110 further includes a receiver 114 for receiving and processing electrical signals received from the transducer, and a transmitter 116 for sending electrical signals to the transducer. The sonar unit 110 further includes a transmitter controller 118 for controlling the operation of the transmitter including the start and stop, and the frequency of a ping.

The signals received by the receiver 114 are sent to a preprocessor for conditioning and compensation. Specifically, the preprocessor 120 includes a filter conditioner 122 for eliminating outlier values and for estimating and compensating for hydrophone variations. The preprocessor further includes a Doppler compensator 124 for estimating and compensating for the motion of the vehicle. The preprocessed signals are sent to a matched filter 130.

The matched filter 130 includes a pulse compressor 132 for performing matched filtering in range, and a beamformer 134 for performing matched filtering in azimuth and thereby perform direction estimation.

The signal corrector 140 includes a grazing angle compensator 142 for adjusting sonar images to compensate for differences in grazing angle. Typically, if a sonar images a collection of point scatterers the image varies with observation angle. For example, a SAS system operating at a fixed altitude and heading observing a sea floor path will produce different images at different ranges. Similarly, SAS images made at a fixed horizontal range would change if altitude were varied. In such cases, changes in the image would be due to changes in the grazing angle. The grazing angle compensator 142 is configured to generate grazing angle invariant images. One such grazing angle compensator is described in U.S. patent application Ser. No. 12/802,454 titled “Apparatus and Method for Grazing Angle Independent Signal Detection,” the contents of which are incorporated herein by reference in their entirety.

The signal corrector 140 includes a phase error corrector 144 for correcting range varying phase errors. The phase error corrector 144 may correct for phase error using a technique described with reference to FIG. 7. Generally, the phase error corrector 144 breaks the image up into smaller pieces, each piece having a substantially constant phase error. Then, the phase error may be estimated and corrected for each of the smaller pieces.

The system 100 further includes a signal detector 150 having a signal correlator 152 and a storage 154. The signal detector 150 may be configured to detect potential targets, estimate the position and velocity of a detected object and perform target or pattern recognition. In one implementation, the storage 154 may include a map store, which may contain one or more previously obtained SAS images real aperture images or any other suitable sonar image. The signal correlator 152 may be configured to compare the received and processed image obtained from the signal corrector 140 with one or more prior images from the map store 154.

The system 100 may include other components, not illustrated, without departing from the scope of the present disclosure. For example, the system 100 may include a data logging and storage engine. In certain implementations the data logging and storage engine may be used to store scientific data which may then be used in post-processing for assisting with navigation. The system 100 may include a security engine for controlling access to and for authorizing the use of one or more features of system 100. The security engine may be configured with suitable encryption protocols and/or security keys and/or dongles for controlling access. For example, the security engine may be used to protect one or more maps stored in the map store 154. Access to one or more maps in the map store 154 may be limited to certain individuals or entities having appropriate licenses, authorizations or clearances. Security engine may selectively allow these individuals or entities access to one or more maps once it has confirmed that these individuals or entities are authorized. The security engine may be configured to control access to other components of system 100 including, but not limited to, navigation controller 170, motor controller 180, sensor controller 190, transmitter controller 118, and CCU 160.

Generally, with the exception of the transducer 112, the various components of system 100 may be implemented in a computer system, such as computer system 200 of FIG. 2. More particularly, FIG. 2 is a functional block diagram of a general purpose computer accessing a network. The holographic navigation systems and methods described in this application may be implemented using the system 200 of FIG. 2.

FIG. 2 is block diagram of an exemplary computer system for implementing at least a portion of the systems and methods described in the present disclosure. The exemplary system 200 includes a processor 202, a memory 208, and an interconnect bus 218. The processor 202 may include a single microprocessor or a plurality of microprocessors for configuring computer system 200 as a multi-processor system. The memory 208 illustratively includes a main memory and a read-only memory. The system 200 also includes the mass storage device 210 having, for example, various disk drives, tape drives, etc. The main memory 208 also includes dynamic random access memory (DRAM) and high-speed cache memory. In operation and use, the main memory 208 stores at least portions of instructions for execution by the processor 202 when processing data (e.g., model of the terrain) stored in main memory 208.

In some implementations, the system 200 may also include one or more input/output interfaces for communications, shown by way of example, as interface 212 for data communications via the network 216. The data interface 212 may be a modem, an Ethernet card or any other suitable data communications device. The data interface 212 may provide a relatively high-speed link to a network 216, such as an intranet, internet, or the Internet, either directly or through another external interface. The communication link to the network 216 may be, for example, any suitable link such as an optical, wired, or wireless (e.g., via satellite or 802.11 Wi-Fi or cellular network) link. In some implementations, communications may occur over an acoustic modem. For instance, for AUVs, communications may occur over such a modem. Alternatively, the system 200 may include a mainframe or other type of host computer system capable of web-based communications via the network 216.

In some implementations, the system 200 also includes suitable input/output ports or may use the Interconnect Bus 218 for interconnection with a local display 204 and user interface 206 (e.g., keyboard, mouse, touchscreen) or the like serving as a local user interface for programming and/or data entry, retrieval, or manipulation purposes. Alternatively, server operations personnel may interact with the system 200 for controlling and/or programming the system from remote terminal devices (not shown in the Figure) via the network 216.

In certain implementations, a system implementing high frequency holographic navigation requires a processor, such as a navigational controller 170, coupled to one or more coherent sensors (e.g., a sonar, radar, optical antenna, etc.) 214. Data corresponding to a model of the terrain and/or data corresponding to a holographic map associated with the model may be stored in the memory 208 or mass storage 210, and may be retrieved by the processor 202. Processor 202 may execute instructions stored in these memory devices to perform any of the methods described in this application, e.g., grazing angle compensation, or high frequency holographic navigation.

The system may include a display 204 for displaying information, a memory 208 (e.g., ROM, RAM, flash, etc.) for storing at least a portion of the aforementioned data, and a mass storage device 210 (e.g., solid-state drive) for storing at least a portion of the aforementioned data. Any set of the aforementioned components may be coupled to a network 216 via an input/output (I/O) interface 212. Each of the aforementioned components may communicate via interconnect bus 218.

In some implementations, a system implementing high frequency holographic navigation requires a processor coupled to one or more coherent sensors (e.g., a sonar, radar, optical antenna, etc.) 214. Examples of suitable sensor arrays are illustrated schematically in FIGS. 3 and 4. An exemplary sonar array is shown in FIG. 3. This array includes a transmitter, receive array, and receive element. An exemplary synthetic aperture sonar array is shown in FIG. 4. This array includes a transmitter, receive array, and receive element, and a virtual array with an associated phase center/virtual element.

Data corresponding to a model of the terrain, data corresponding to a holographic map associated with the model, and a process for grazing angle compensation may be performed by a processor 202 operating on the data, as shown in FIG. 2. The system may include a display 204 for displaying information, a memory 208 (e.g., ROM, RAM, flash, etc.) for storing at least a portion of the aforementioned data, and a mass storage device 210 (e.g., solid-state drive) for storing at least a portion of the aforementioned data. Any set of the aforementioned components may be coupled to a network 216 via an input/output (I/O) interface 212. Each of the aforementioned components may communicate via interconnect bus 218.

In operation, a processor 202 receives a position estimate for the sensor(s) 214, a waveform or image from the sensor(s) 214, and data corresponding to a model of the terrain, e.g., the sea floor. In some implementations, such a position estimate may not be received and the process performed by processor 202 continues without this information. Optionally, the processor 202 may receive navigational information and/or altitude information, and a processor 202 may perform a coherent image rotation algorithm. The output from the system processor 202 includes the position to which the vehicle needs to move to.

The components contained in the system 200 are those typically found in general purpose computer systems used as servers, workstations, personal computers, network terminals, portable devices, and the like. In fact, these components are intended to represent a broad category of such computer components that are well known in the art.

It will be apparent to those of ordinary skill in the art that methods involved in the systems and methods of the invention may be embodied in a computer program product that includes a non-transitory computer usable and/or readable medium. For example, such a computer usable medium may consist of a read only memory device, such as a CD ROM disk, conventional ROM devices, or a random access memory, a hard drive device or a computer diskette, a flash memory, a DVD, or any like digital memory medium, having a computer readable program code stored thereon.

Optionally, the system may include an inertial navigation system, a Doppler sensor, an altimeter, a gimbling system to fixate the sensor on a populated portion of a holographic map, a global positioning system (GPS), a long baseline (LBL) navigation system, an ultrashort baseline (USBL) navigation, or any other suitable navigation system.

High-Frequency Holographic Navigation

FIGS. 5A-5B depict processes 500 and 550 for navigating a terrain using an exemplary high-frequency sonar navigation system, such as system 100. In particular, the processes 500 and 550 may be implemented across several components of system 100 of FIG. 1. The system 100 may receive via wire or wirelessly, at the map store 154, a prior high frequency SAS image of a portion of the terrain being navigated (step 502). The prior image may have been obtained using a frequency range greater than 100 kHz. For example, the frequency of the prior SAS image may include a well-formed image in the frequency range of 100 kHz-110 kHz, or between 110 kHz-120 kHz. The frequency of the prior SAS image may be between 100 kHz-150 kHz and/or 125 kHz-175 kHz and/or 175 kHz-225 kHz. The frequency of the prior SAS image may be greater than 500 kHz and in certain implementations, the frequency may be greater than 1 MHz. In certain implementations, the frequency ranges may be selected based on application. For example, for certain ocean systems, the frequencies may be up to about 500 kHz, and in certain medical ultrasound systems, the frequency may be about 15 MHz. In certain implementations, the frequency ranges may be selected to be less than 100 kHz. In such implementations, the process 500 may be especially beneficial depending on the ratio of the size of the error to the wavelength. In one example, for ships bouncing in waves process 500 may be beneficial for frequencies down below 10 kHz. The prior SAS image may be grazing angle compensated and/or phase error corrected and the frequency of the image may be post-grazing angle correction. In certain alternative implementations, the prior image may include a low frequency image in the tens of kHz and less than 100 kHz.

The process 500 includes predicting, by the CCU 160, an initial position value of the vehicle traversing the terrain based on a previous position (Step 504). The position value may be represented in any suitable coordinate system. The CCU 160 may generate this initial position value based on information from the navigational controller 170 and previous motion. The CCU 160 may also determine an error estimate or navigational uncertainty associated with this initial position value (step 506).

The system 100 may insonify a portion of the terrain being navigated with a high-frequency signal and generate a current sonar image (step 508 in FIG. 5A, step 558 in FIG. 5B). In certain implementations, CCU 160 in connection with the transmitter controller 118 may send transmission instructions to the transmitter 116 and the transducer array 112. To allow for coherent correlation, the frequency of imaging may be selected to overlap with the frequency range of the received prior map obtained in step 502. The overlap in frequencies may be a complete overlap, a partial overall or an implicit overlap. In a complete overlap, the frequency range of the current sonar image may lie completely within the frequency range of the prior map. In a partial overlap, the frequency range of the current sonar image may partially overlap with the frequency range of the prior map. Even when frequency ranges of the current raw imaging process and the prior mapping process do not overlap, there may still be an implicit overlap if the aspect or viewing angle of the two images are appropriately different. In such an implicit overlap scenario, the grazing angle compensated frequencies of the current image and the prior map at least partially overlap. For example, a 100 kHz signal at a 45 degree grazing angle would have the same projected wavelength as a 70.7 kHz signal at a grazing angle of zero and would consequently constitute an implicit overlap. As another example, a 100 kHz-110 kHz image at a relatively flat grazing angle may coherently correlate with a 110 kHz-120 kHz map at a relatively steeper grazing angle. Generally, the system 100 may operate as a SAS and obtain a high-frequency SAS image. To obtain a SAS image, the system 100 may operate in sectorscan mode, sidescan mode, stripmap mode and/or spotlight mode. System 100 may even operate as a phased array and obtain a real aperture image of the terrain.

The obtained current image, which may include a real aperture image, in step 508 in FIG. 5A or a synthetic aperture image, in step 558 in FIG. 5B, may be passed through preprocessor 120, matched filter 130 and received at the signal corrector 140. The obtained current image is modified to compensate for grazing angle by the grazing angle compensator 142 (step 509 in FIG. 5A, step 559 in FIG. 5B). Generally, the obtained current image is converted to a grazing angle invariant image. The grazing angle compensator 142 approximates the terrain (e.g., sea floor) as a smoothly undulating manifold with embedded point scatterers, and models the sonar signals as interference between point scatterer echoes. Shadowing and occlusion are generally neglected and changes in grazing angle are assumed to change the pitch of the echo. Changes in pitch generally cause all frequencies to be scaled by a multiplier which is the secant of the grazing angle. By reversing the process (i.e., projecting the echo onto the sea floor), a relationship between scatterer spacing and image frequency is established that is independent of grazing angle. Typically, grazing angle compensation is limited by transmitter design; the applicable range of angles is determined by signal bandwidth and transmitter properties.

The compensated obtained current image, which may be a real aperture image or a synthetic aperture image or any suitable sonar image, is then modified to correct for range varying phase errors by the phase error corrector 144 (step 510 in FIG. 5A, step 560 in FIG. 5B). The process of correcting for range varying phase errors, which allows for high-frequency imaging and navigation, is described in more detail with reference to FIG. 7.

The compensated and error corrected obtained current image of the terrain is coherently correlated, at the signal correlator 152, with the prior SAS map received at step 502 (step 511 in FIG. 5A, step 561 in FIG. 5B). Generally, because image intensities can spatially vary, the signal detector 150 may be configured to perform a normalized correlation. In certain implementations, the normalized correlation may be performed by calculating the correlation coefficient. Generally for sonar images, the correlation coefficient is often low (less than 0.1) and the values depend on the available structure. Without a priori knowledge of the terrain, it is difficult to define detector thresholds. Detection may still be possible, however, because signal to noise ratios (SNR) may be high. The signal detector 150 may calculate additional statistics of the normalized correlation include the statistical distributions of the signals (amplitude and phase) and/or noise. The statistical distributions may include Rayleigh and/or Gaussian distributions. The detector thresholds may be selected based on the distribution. Examples of suitable correlation techniques included in signal detector 150 techniques described in “On Correlating Sonar Images,” Richard J. Rikoski and J. Tory Cobb and Daniel C. Brown, Robotics: Science and Systems '05, 2005, and “Holographic navigation,” Richard J. Rikoski and Daniel C. Brown, ICRA '08, 2008, the entire contents of each of which are incorporated herein by reference.

Using the coherent correlation of the image with the map, the CCU 160 may determine a measured position value (step 512) and the associated error estimate of the position (step 513). Certain exemplary techniques to determine measured position are described in “Holographic navigation,” Richard J. Rikoski and Daniel C. Brown, ICRA '08, 2008, the entire contents of which are incorporated herein by reference.

Based on the new position estimate, the CCU 160 may update control signals and instruct the motor controller 180 to move the vehicle accordingly. If navigation to a new location is complete (step 515), the previous map-based calculated position is set as the previous position (step 516) and the process 500 is repeated at the new position.

In some implementations, traditional SAS navigation may robustly and easily solve for position, but may less efficiently solves for heading, e.g., of an unmanned autonomous vehicle. It may be possible to correlate a synthetic aperture image against either a real aperture or synthetic aperture image at various angles to estimate the heading, but this may be computationally intensive. The proposed system may solve this problem by decomposing the heading estimation problem into a two step process. First, holographic navigation is used to estimate position. Then, a correlation is performed using an angular coordinate system centered on the estimated position. Assuming range is r and angle is θ, a coordinate system which is a function of r and θ is used for correlation. In the simplest instantiation, (f(r), g(θ))=(r,θ), but alternatives like f(r)=horizontal distance along the bottom or g(θ)=sin(θ) may also be appropriate. The correlation may either be in range and angle or just in angle, but to detect heading it may be necessary to correlate in angle. The angle with the maximum correlation corresponds to the direction the reacquisition sonar is facing. Two exemplary processes for solving for heading are processes 600 and 650 illustrated in FIG. 6A and FIG. 6b , respectively. Processes 600 and 650 may be similar to processes 500 and 550, respectively, except for the step of determining heading based on measured position value as shown in step 602 in FIG. 6A and step 652 in FIG. 6B.

As noted above, the system 100 includes a phase error corrector 144 to correct for range varying phase errors. Range varying phase errors may lead to low correlation values when correlating between two high frequency images (e.g., an image and a map). As an illustrative example, suppose a robot with a high frequency sonar attempts to correlate its imagery with a prior map but it has a 1 cm altitude error and a 1 cm wavelength. Directly underneath the vehicle this leads to a 2 cm path length error, or 2 cycles. At long range, this leads to a zero cycle delay. Consequently, if the prior sonar image and the conjugate of the new sonar are multiplied together but not summed (image1*conj(image2)) what will be observed is a range varying phase that is due to that altitude error. When the multiplied images are summed together this range varying phase error will cause destructive interference, leading to a very low correlation value.

Similarly, as a second illustrative example, assume a 0.01% sound speed error, a 1 cm sonar, and an operating altitude of 5 meters. The travel path directly under the vehicle is 10 meters, or 1000 cycles, leading to a 1/10^(th) of a cycle error. At a range of 50 meters (or 10,000 cycles) this leads to a full cycle error. At a range of 500 meters this leads to 10 cycles of error. So again, when the multiplied sonar images are summed constructive interference will drive the cross correlation down. High frequency holographic navigation attempts to solve these problems by either using image pieces which are small enough to be immune to those effects, or by estimating and correcting for those biased parameters.

One method of correcting for phase error is to cut the image into small regions where the phase error is constant and use those as independent measurements. In regions with very high signal-to-noise ratio (SNR) this may be very efficient.

FIG. 7 depicts a process 700 for correcting range varying phase errors in a high-frequency sonar system. The process 700 may be implemented on phase error corrector 144 of the signal corrector 140 in system 100 of FIG. 1. The process 700 begins with receiving a real aperture image or synthetic aperture image.(step 702). In certain implementations, the real aperture image may be modified with grazing angle compensation. The phase error corrector 144 may estimate the range varying phase error of the entire received real aperture image (step 704). The phase error corrector 144 may then determine if the variation in phase error across the image is less than an error threshold (step 706). The error threshold may be set as desired. In certain implementations, the error threshold may be set depending on the maximum range of the real aperture image. If the variation in phase error across the image is greater than the error threshold then the image may be split into sub-regions (step 708). The phase error corrector 144 may split the image into sub regions as desired. The sub-regions may be of equal sizes or of different sizes. Sub-regions may be of varying sizes such that the size variation may be based on the range. The phase error corrector 144 may estimate the phase error for each sub-region (step 712), and determine if the variation in phase error across each subregion is less than a phase error threshold (step 714). The error threshold for subregions may be the same as or different from (greater than or less than) the error threshold associated with step 706. The error thresholds for each subregion may be different or the same. If the phase error in a particular subregion is less than the error threshold, that particular subregion may be corrected for the corresponding phase error, which is substantially constant across the entirety of the particular subregion (step 715). If the phase error is greater than the threshold then the subregion may be split into smaller subregions and steps 708, 712, 714 may be repeated. In certain implementations, the subregions may be selected such that they have constant altitude phase error or constant sound speed error. One or more selected subregions in one or more sensitive mapping regions may be selected to have the largest possible size.

In some implementations, suppose system 100 can correlate small patches (e.g., 50 pixels by 50 pixels). In such example, if the image is 1000 by 1000 pixels, then system 100 may cut the image up into 20×20 regions of 50×50 pixels each. The system 100 may perform 400 separate correlations. Each correlation may have a peak with a slight shift and a slightly different phase value due to the unknown error function. System 100 may take the absolute value of each correlation and sum them all together to eliminate destructive interference due to the phase differences. Such an approach may be advantageous at least when the error function is unknown. Also, although the noise may be Rayleigh distributed when the distribution of the absolute value for a single image correlation is viewed, but when system 100 sums a large number together the law of large numbers applies and the noise becomes Gaussian distributed.

The system 100 may include other methods for compensating for phase errors. In certain implementations, the real image is taken and a sum of the envelopes of small image correlation regions with approximately stationary phase is calculated before calculating a probability density function based on that sum. This may be similar to a speckle reduction technique used in imaging methods. The sums can either be for a single altitude solution or for multiple altitude solutions; if multiple altitude solutions are used then the technique measures altitude bias. In certain implementations, using the envelope only (the absolute value of the correlation result) removes the relative phase differences between correlation results. It is important to note that summing together a large number of correlation images results in a transition from Rayleigh to Gaussian distributed speckle intensity; this difference may be important when converting the correlation result to a probability density function. Previous holographic navigation techniques, which use the Rayleigh distribution, typically fail when presented with correlation results based on sub-image summation; switching to a more representative distribution is key. When a small number of images are summed together the distribution may not yet be fully Gaussian and may be better represented by some other distribution such as a K-distribution.

In some implementations, patches may be used with roughly stationary phase to estimate the range varying phase error and then apply an appropriate correction so as to enable full waveform correlation. In some implementations, estimating the range varying phase error may be done several ways, including, inter alia, unwrapping the phase and fitting a curve, doing a least squares fit to the raw angles, or changing coordinate systems and using a Fast Fourier Transform (FFT) or any type of fourier transform such as a Discrete Fourier Transform (DFT) or a wavelet transform to find the delay. This implementation (combining a coordinate system change and an FFT) is applicable to time delay estimation beyond holographic navigation (for instance, to motion estimation for synthetic aperture systems using displaced phased center navigation, especially for heave estimation).

In one example of a heave estimation technique includes the following:

$t = {\frac{2r}{c} = \frac{2\sqrt{x^{2} + z^{2}}}{c}}$ Where t=time, r=range, c=sound speed, x=horizontal range and z=altitude.

$x = \sqrt{\frac{c^{2}t^{2}}{4} - z^{2}}$ $\frac{dx}{dz} = {{- \frac{z}{\sqrt{\frac{c^{2}t^{2}}{4} - z^{2}}}} = {- \frac{z}{x}}}$ ${x^{\prime} \cong {x + e_{x} + {\frac{dx}{dz}e_{z}}}} = {x + e_{x} - {\frac{z}{x}e_{z}}}$ Where e_(x)=x error, e_(z)=z error, s˜=envelope function, k_(x0)=wavenumber of carrier frequency.

Convert s(t)−>s(x)

Now, s1(x1)!=s2(x2) because of error (e_(x), e_(z)). So when you multiply s1(x1) by the conjugate of s2(x2), the signals differ by an offset:

${x\; 2} = {{x\; 1} + e_{x} - {\frac{z}{x\; 1}e_{z}}}$ ${{x\; 2} - {x\; 1}} = {e_{x} - {\frac{z}{x\; 1}e_{z}}}$

Now, e_(x) is usually pretty easily observable and removable, getting you to a range varying error:

${{x\; 2} - {x\; 1}} = {{- \frac{z}{x\; 1}}e_{z}}$

The signals are now

s 1(x) = (x)e^(jk_(x 0)x) ${s\; 2(x)} = {{(x)e^{{jk}_{x\; 0}x}} = {\left( {x + {\frac{dx}{dz}e_{z}}} \right)e^{{jk}_{x\; 0}{({x + {\frac{dx}{dz}e_{z}}})}}}}$ ${s\; 1(x)s\; 2^{*}(x)} = {(x)(x)e^{{- {jk}_{x\; 0}}\frac{dx}{dz}e_{z}}}$ ${s\; 1(x)s\; 2^{*}(x)} = {(x)(x)e^{{jk}_{x\; 0}\frac{z}{x}e_{z}}}$

Now, changing into a new coordinate system

$\sigma = \frac{k_{x\; 0}z}{x}$ changes this into s1(σ)s2*(σ)=

(σ)

(σ)e ^(jσe) ^(z) ≅

(σ)² e ^(jσe) ^(z) From here, we can do a Fourier transform to estimate e_(z).

In certain implementations, one or more techniques employed by the signal corrector 140 or any other component of system 100 may be used in certain applications including, but not limited to, change detection and two pass interferometry. Change detection is typically a process of taking two passes by a scene, then accurately aligning two images, and coherently comparing them. In regions where there has been “change” they may decorrelate significantly. Two-pass interferometry is typically a process of taking two passes over a scene, aligning two images and comparing the phase of these two images. The comparison of phase of two images may reveal changes and deformation in the terrain. In certain implementations, such deformation may be over timespans of days to years. Such applications may be useful for geophysical monitoring of natural hazards, such as earthquakes, volcanoes and landslides, and also in structural engineering, including monitoring of subsidence and structural stability. Other applications of system 100, and particularly signal corrector 140, include reconnaissance, surveillance and targeting. These applications may use system 100 to generate high resolution images and to distinguish terrain features (surface and/or underwater) and to recognize and identify selected man made targets. Still other applications include interferometry, navigation, guidance, imaging foliage and underground or subsurface targets, detecting and moving targets, and environmental monitoring application such as monitoring oil spills.

They systems and methods described herein may be adapted as desired for both sonar and radar systems. For example, sonar transducers may be replaced with suitable radar transducers, and one or more components may be modified, added to or removed from the systems described herein to operate in a sonar and radar regime. In some implementations, the systems and methods may be configured to operate as both sonar and radar devices, without departing from the scope of the present disclosure. In certain implementations, when the systems and methods are configured for sonar imaging, the frequencies may be in the range from 100 kHz to about 200 kHz. In certain implementations, when the systems and methods are configured for radar imaging, the frequencies may be in the range from 1 GHz to about 30 GHz. Generally, the systems and methods described herein may be applied for any frequency range, without departing from the scope of the present disclosure.

Certain Applications of High-Frequency Holographic Mapping and Navigation Systems

The systems and methods described here make use of various other aspects of the holographic nature and high frequency of synthetic aperture images, which inventor has recognized. For example, systems and methods are described herein for determining a three-dimensional model of a shape based on its two dimensional shading and shadowing of acoustic signals. In some implementations, traditional shape from shading may be one parameter shy of a solution; an approximation may be required in order to derive a three dimensional model from a pure image. However, since a synthetic aperture image is a quasi-hologram and contains a continuum of images of a range of angles, it may contain enough information to over-constrain the shape from shading problem. The systems and methods above described solve the shape from shading problem by decomposing the SAS image into lower resolution sub-patches and then deriving their orientation from the shading observed from multiple vantage points.

The systems and methods described herein include methods for positioning sensors (such as Tsunami sensors) and navigation beacons with high-precision using HF holographic navigation. In certain implementations, Tsunami buoys use sensors on the seafloor to small variations in water pressure. To make accurate measurements it is necessary that the sensor be positioned properly on the seafloor. If the sensors are hidden behind rocks or are not level/well placed on the sea floor it can affect the accuracy of their measurements. The systems and methods described herein combine a holographic navigation system, and maneuvering system, and a tsunami sensor so that tsunami sensors can be positioned very precisely using a prior map.

In some implementations, holographic navigation enables very high precision navigation relative to the seafloor, but may be limited by the necessity of periodically observing the seafloor. Midwater system that may be at a great altitude to observe the seafloor cannot take advantage of holographic navigation or its precision. The systems and methods described herein address this limitation by combining a beacon system, and maneuvering system, and a holographic navigation system. The beacon is able to position itself very precisely, enabling systems including long baseline navigation or ultrashort baseline navigation without needing to calibrate the beacon system using a ship.

The systems and methods described herein include methods for monitoring and modeling a water column using an autonomous underwater vehicle (AUV) based on high-precision location measurements obtained using HF holographic navigation. In some implementations, AUVs either circle a buoy or simply form a wagon wheel. By transmitting orthogonal signals to one another they can measure time of flight between positions and also measure differential time of flight. From time of flight, it may be possible to determine the sound speed of the water; from differential time of flight it may be possible to determine the Doppler shift along the connecting vector/estimate water velocity. Vehicle positions are determined using holographic navigation, thereby enabling a high precision model of the water column in post processing. Vehicles may dock at a central buoy for recharging.

In some implementations, AUV vehicle recovery and vehicle docking may be difficult problems due to the dynamic nature of both the vehicle and the destination. If it can be decomposed into a purely relative problem, the vehicle needs its position relative to the dock as well as its orientation. The systems and methods described herein allow a vehicle to passively estimate its non-range position relative to a docking system passively, and allows it to estimate its range to the docking system actively.

In some implementations, the system takes advantage of the fact that a blazed array transmits different frequencies at different angles. Using two blazed arrays with different frequencies oriented orthogonally it creates a two dimensional grid of frequencies. For instance, suppose a 300-600 kHz blazed array was oriented such that the frequencies varied with horizontal displacement, and a 600-1200 kHz blazed array was oriented such that its frequencies varied in the vertical direction. A vehicle observing 450 kHz and 900 kHz would be driving straight into the dock. A system observing 500 kHz and 900 kHz would have the correct elevation but would be displaced horizontally. A system observing 450 kHz and 950 kHz would be displaced vertically. In some implementations, to measure the vehicle's orientation with respect to the docking station, the vehicle would have a small passive array to measure the direction of the incoming signal from the blazed arrays.

In some configurations, range may be measured using a small beacon system such as an ultra short baseline beacon. In some implementations, range may be measured using high frequencies that are only observable at short ranges, or may be neglected entirely (purely a glide path based docking method).

Seismic Survey System Using Planar SAS and Holographic Navigation

Seismic survey is generally a form of 2D or 3D geophysical survey that is used to measure terrestrial or extra-terrestrial properties by means of acoustics or electromagnetic. Seismic survey systems are necessary for offshore oil exploration, but they are large, ship intensive, expensive, and high power. Traditional seismic survey systems use very high powered transmitters to insonify the bottom, and receive the signal on a network of towed arrays which are dragged behind a large ship.

In certain implementations, the systems and methods include a seismic survey system having a combination of orthogonal transmitters and multiple receivers to form a full planar synthetic aperture sonar with higher resolution, lower power, and fewer large ships than a traditional seismic survey system.

Applicants' system takes advantage of the phase center approximation of synthetic aperture sonar (SAS). A phase center is located halfway between the transmitter and receiver. For an array to be fully populated (from a Nyquist perspective) it needs to have an appropriate number of properly spaced phase centers.

In certain implementations, system 100 includes multiple transmitters. Using multiple transmitters with orthogonal signals, it may be possible to distinguish between phases created by different transmitters. Therefore, by using M transmitters and N receivers, it is possible to create MN phase centers. This is often less expensive than using one transmitter and MN receivers. In certain implementations, the system 100 generalizes to any practical value of M and N.

The transmitters of the system can be mounted on any sort of vessel or robot (ship, autonomous underwater vehicle (AUV), unmanned surface vehicle (USV), nuclear submarine, etc). In certain implementations, the transmitters of the systems described herein may require a relatively high power. In such implementations, the vessel may be equipped with suitable power delivery systems to supply the needed power to the transmitters. One example of a vessel includes modest sized USVs such as 10 m RHIBs (Rigid Hull Inflatable Boats) since autonomous systems are ideal for maneuvering in formation and surface craft enable the use of GPS.

In certain implementations, system 100 includes multiple receivers. The receivers of the system can be mounted on various vehicles and use various array types without dragging arrays behind ships (even though this is possible). In some implementations, AUVs are flown in formation close to the bottom. This reduces losses on the return path, reduces the necessary transmit power, and allows the receivers to be precisely positioned using holographic navigation. Holographic navigation in this manner requires a prior seabed survey of the area, but this is relatively inexpensive.

In some implementations, the combination of transmitters and receivers form a line array of phase centers. That line array is then translated orthogonally to its axis in a predominantly horizontal direction, so that when the data is accumulated there is a Nyquist sampled planar array. Using that planar array it may be possible to beamform the signal to form a 3D image composed of high resolution voxels penetrating deep into the seabed.

Generally, for seismic survey applications, system 100 may be operated at any suitable frequency without departing from the scope of the present disclosure. For example, system 100 may be configured for frequencies in the range of 1 Hz to about 10 kHz. The system 100 may be generally configured for frequencies less than 10 kHz, including frequencies in the range of 100 Hz to about 10 kHz. In certain implementations, system 100 may be adapted with electromagnetic transducers and suitable components for radar-based seismic applications. In such applications the frequencies may range from about 300 MHz to about 30 GHz. Whether configured to operate for radar or sonar based applications, system 100 may use frequencies in any suitable range without departing from the scope of the invention.

They systems and methods described herein may be adapted as desired for both sonar and radar systems, and accordingly for both synthetic aperture sonar (SAS) and synthetic aperture radar (SAR) systems. For example, sonar transducers may be replaced with suitable radar transducers, and one or more components may be modified, added to or removed from the systems described herein to operate in a sonar and radar regime. In some implementations, the systems and methods may be configured to operate as both sonar and radar devices, without departing from the scope of the present disclosure. In certain implementations, when the systems and methods are configured for sonar imaging, the frequencies may be in both high and low frequency ranges in the range from 10 Hz to about 1 kHz. In certain implementations, when the systems and methods are configured for radar imaging, the frequencies may be in the range from 1 MHz to about 100 MHz. Generally, the systems and methods described herein may be applied for any frequency range, without departing from the scope of the present disclosure.

Low Grating Sidelobe SAS

The systems and methods described herein include adding multiple transmitters and generating orthogonal pinging sequences configured to enhance the performance of a SAS system. In particular, the systems and methods described herein include a SAS having a low-grating sidelobe (as described with reference to FIGS. 8-9B), a SAS having a high coverage rate using multiple transmitters (as described with reference to FIGS. 10-11C), and an overpinging sequence for increasing the range of the SAS system.

In general, grating sidelobes occur when active sonar elements are one or more wavelengths apart. Grating sidelobes may not be fully suppressed when elements are more than a half wavelength apart. For most active sonar systems this spacing is impractical, since it would require an extremely high channel count and omni-directional elements. Instead, most systems use larger transducer elements with limited beam patterns. The resulting beam pattern of the sonar system is the product of the array beam pattern (including grating lobes) and the beam pattern of the individual elements. Since those elements output relatively little in the direction in the direction of the grating lobes, this spacing partially suppresses the lobes. According to one illustrative implementations, the transmitter and receiver elements have a relevant dimension (e.g, width) d, and assuming a phase center approximation (the ellipsoidal travel path between transmitter and receiver is modeled as a range circle centered halfway between the transmitter and receiver), the classical SAS array has “d/2” spacing, as shown in FIG. 9A.

According to the illustrative implementation of FIG. 9B, one may arbitrarily increase the array sampling from d/2 to something higher (e.g., d/4), as shown. This spacing takes advantage of the fact that a synthetic aperture sonar system is in constant motion so that it can accumulate many pings/phase centers so as to create a very high resolution image. The system then operates by transmitting orthogonal signals after a very short delay to create additional phase centers in between the original phase centers so that they can be added in processing to create an array with d/2N sampling, where N>1. In effect, the SNR is no longer grating lobe limited.

FIG. 8 depicts a process 800 for using a plurality of orthogonal signals in a synthetic aperture sonar (SAS) system, such as system 100, to generate images. In particular, process 800 may be configured to use orthogonal signals to generate SAS beams having suppressed grating lobes. Process 800 begins with providing a SAS array (such as array 112 of FIG. 1) having a transmitting element and plurality of receiving element (step 802). Such an array is depicted in FIGS. 9A and 9B. In certain implementations, each of the transmitting and receiving elements may have a first dimension, d. The dimension may include any suitable dimension including length, width and diameter. A user or CCU 160 may determine the number of orthogonal pings, N (step 804). As noted above generally N>1. In certain implementations, N=2, such that the sampling is about d/4. Each ping p(i)={p1, p2 . . . , pn} has a duration of Tp. In certain implementations, each ping may be orthogonal to one or more previous pings such that pings overlapping in time are orthogonal to each other and non-orthogonal pings do not overlap with each other. The CCU 160 or the transmitter controller 118 may calculate a first effective spacing D=d/2, representative of an effective distance between transmitter and receiver elements during motion.

According to process 800, the transmitter controller 118 instructs the transmitter 116 to set time to t0 and start transmitting the ping p(i), where i=1 (step 808). The transducer array 112 is moved along an axis parallel to that connecting the receiving elements (step 810). In certain implementations, it may be acceptable for the face of the transmitting elements or projectors to not be coplanar to the face of the receivers. For example, streamlined vehicles include a polyurethane coating continuous with the body form, however the actual transmitting elements may be embedded about 1-2 inches behind that polyurethane window. The CCU 160 determines if the transducer array 112 has moved a distance of D/N (step 812). If the transducer array 112 has moved a distance of D/N, then the CCU 160 determines if all the pings have been transmitted in the current iteration (step 814). If all the pings have not been transmitted then, the next ping (which is orthogonal to the previous ping) is transmitted and the process is repeated from step 810. For example, if N=2, the transmitter width is d=0.1 m, the robot is translating at 1 meter per second, the time period between the first ping in each cycle is 1 s, and the pulse length is 0.2 s, one implementation of the process 800 includes firing the first ping at time t=0. At time, t=0.025 s, the transducer 112 may have translated a distance of d/4. At time, t=0.025 s, the transducer 112 may be configured to fire the second orthogonal ping. Between 0.025 s and 0.2 s, both the first ping and the orthogonal second ping are transmitting. Between 0 s and 0.025 s, only the first ping is transmitting, and between 0.2 s and 0.225 s, only the second orthogonal ping is transmitting.

Generally, and not to be bound by theory, the process allows for delaying the transmission of the second signal until the vehicle has translated enough to create a second virtual array. As was shown in process 800, the delay may be related to vehicle speed and firing may be adjusted based on the measured motion while keeping the vehicle speed constant (“slave to speed” configuration). In certain other implementations, the delay may be fixed and the vehicle speed may be adjusted, including performing alongtrack compensation. In another configuration, the matched filter length may be adjusted slightly to compensate for alongtrack motion imperfections when defining phase centers (e.g. in the above example transmitting a noise sequence, but then dropping the first 0.001 s to 0.00001 s of the matched filter template to correctly place the effective vehicle ping start position.

Consider an exemplary system with a 1 meter long broadside array consisting of 10 centimeter elements and a 10 cm transmitter. The virtual array of phase centers is then 50 cm long with phase centers spaced 5 cm apart. In a typical SAS, the vehicle would transmit, move 50 cm, and transmit again. For a variant of the present disclosure operating with d/4 spacing, the vehicle would transmit, move 2.5 cm, transmit an orthogonal signal (so as not to jam the original signal), move an additional 47.5 cm, and then repeat. If the vehicle was moving at 1 m/s, the delay between signals would be 1/40 of a second, or 25 ms. If the transmit signal is longer than 25 ms then the two orthogonal signals will overlap. In this case, the signals need to be designed such that when they are summed together they do not saturate the transmitter. If the goal is higher sampling that d/4, it may be necessary to sum together multiple signals.

This method is not only restricted to broadside synthetic aperture sonar. Broadside active phased arrays may use this technique to form very short aperture synthesis to reduce sidelobes (i.e. a sidescan sonar would fuse two pings). Real aperture and synthetic aperture forward looking and/or squinted sonars may use the same technique to increase their element count. The technique would work very well with circular SAS arrays.

Then the signal can be changed from ping to ping to further reduce sidelobes after aperture synthesis and to suppress noise from the orthogonal signals. This is manifested in several ways, including. Changing a ping changes its autocorrelation function during aperture synthesis; summing together different autocorrelation functions with different sidelobe structures will reduce the relative magnitude of those sidelobes. Changing the ping changes the cross correlation function between the two subcomponent pings so that during aperture synthesis, the noise is not locally a standing wave and instead destructively interferes.

They systems and methods described herein may be adapted as desired for both sonar and radar systems. For example, sonar transducers may be replaced with suitable radar transducers, and one or more components may be modified, added to or removed from the systems described herein to operate in a sonar and radar regime. In some implementations, the systems and methods may be configured to operate as both sonar and radar devices, without departing from the scope of the present disclosure. In certain implementations, when the systems and methods are configured for sonar imaging, the frequencies may be in both high and low frequency ranges in the range from 10 kHz to about 200 kHz. In certain implementations, when the systems and methods are configured for radar imaging, the frequencies may be in the range from 100 MHz to about 30 GHz. Generally, the systems and methods described herein may be applied for any frequency range, without departing from the scope of the present disclosure.

High Coverage Rate SAS

In certain implemenations, the present disclosure relates to a device for a synthetic aperture sonar with a real array with N elements of size d which combine to create a real array of length L. Because of the phase center approximation, the effective position of the elements is halfway between the transmitter and receivers, making the effective array length L/2. This effective array will be referred to as a virtual array, as shown in FIGS. 11A-C. In some implemenations, if two vertically displaced transmitters are used which transmit orthogonal signals, it may be possible to create two vertically displaced virtual array and perform interferometry. One example of vertically displaced transmitters is described in U.S. Pat. No. 8,159,387, entitled “Multi-transmitter Interferometry,” the contents of which are incorporated herein by reference in their entirety.

In some implementations, if M vertically displaced transmitters are used, it is possible to create M virtual arrays. In some implementations, if two horizontally separated transmitters spaced L apart are used, it is possible to create two abutting virtual arrays giving the vehicle an effective array length of L. One example of a multi-transmitter array is described in U.S. Pat. No. 5,295,188 entitled “Synthetic Aperture Side-Looking Sonar Apparatus,” the contents of which are incorporated herein by reference in their entirety. In some instances, SAS may use two transmitters placed away from the receive array to achieve this effective array length. Doubling the array length is generally desirable since it doubles the area coverage rate of a SAS. (Since the vehicle moves one effective array length between pings, if the array length doubles the range of the sonar doubles. If the robot maintains its ping rate, it must double its velocity in order to be in position for the subsequent ping. In either case, the coverage rate doubles.)

Inventor has recognized that placing a pair of transmitters away from the array spaced L apart results in an effective sonar array length of L (which is greater than prior art systems that have an effective array length of L/2). The transmitters are typically placed away from the array because if they are placed on either side on the receive array they will not be L apart, but L+D. This spacing results in virtual arrays which may have a missing element, resulting in grating sidelobes.

Inventor's method allows for a more flexible placing of transmitters, allowing for a larger area coverage rate in a smaller package. The method includes using multiple transmitters with orthogonal signals fired non-synchronously and using delays and vehicle translation to form abutting virtual arrays. For a simple array of N elements of size d, with transmitters of size d on either side of the array, the forward transmitter will start transmitting first, followed by a delay as the vehicle moves d/2 forward, then the aft transmitter starts firing as shown in FIGS. 11A-C. If, for technical reasons, the transmitters have a different spacing, the timing may be adjusted accordingly.

In some implementations, the method allows for a multitude of transmitters placed along the vehicle. For instance, if four transmitters were used spaced L apart (for a total length of 3L), the effective array length is 2L, and the area coverage rate of the system quadruples over a baseline SAS. In some implementations, if M transmitters are used to lengthen the array, the area coverage increases to M times the baseline coverage. Likewise, pairs (or larger sets) of vertically displaced transmitters may be added to create a longer interferometric array. In the event that the separated along track transmitters cannot be placed in the same vertical position, resulting in virtual arrays that are parallel but not collinear, grazing angle compensation can be used to correct for the vertical displacement. Changing the orthogonal signals used by the transmitters from ping to ping may further reduce noise suppression, which may be desirable. In certain implementations, it may be desirable to fire channels nearly simultaneously to lengthen the array. In such implementations, projectors may be positioned at 0 (e.g., the fore end of the receiver array) and L+w (e.g., the aft end of the receiver array), and other places at approximately (k+/−delta)*L where k is an integer and delta is some acceptable variation. Delta may correspond to a delay that is less than the time before the bottom bounce.

FIG. 10 depicts a process 1000 for transmitting pulses from a SAS system having such a multiple transmitter arrangement. The process 1000 begins with proving a transducer array having a receiver array with a plurality of receiver elements and two transmitter elements (step 1002), each having a width, w. In certain implementations, the transmitting elements are positioned on either side of the receiver array and along the axis of the receiver array, such as the array depicted in FIGS. 11A-C. Transmitter T1 may be position in the aft position and transmitter T2 may be positioned in the fore position such that the vehicle moves in a fore direction. The transmitter T2 may ping first (step 1004), after which the vehicle may move fore alongtrack and along the array axis connecting T1 and T2 (step 1006). The sonar system 100 may determine if the sonar array has moved a distance w/2 (step 1008). Once the sonar system 100 has moved such a distance, the transmitter T1 may ping, such that T1's ping is orthogonal to T2's previous ping (step 1010). The vehicle may then move fore alongtrack (step 1012), and the sonar system 100 may once again query whether the distance traveled is equal to w/2 (step 1014). If such a distance has been traversed, the process 1000 may be repeated from step 1004 and T2 may ping again. If the transmitters of width w were spaced further from the array then the delay would need to be slightly greater.

They systems and methods described herein may be adapted as desired for both sonar and radar systems, and accordingly for both synthetic aperture sonar (SAS) and synthetic aperture radar (SAR) systems. For example, sonar transducers may be replaced with suitable radar transducers, and one or more components may be modified, added to or removed from the systems described herein to operate in a sonar and radar regime. In some implementations, the systems and methods may be configured to operate as both sonar and radar devices, without departing from the scope of the present disclosure. In certain implementations, when the systems and methods are configured for sonar imaging, the frequencies may be in both high and low frequency ranges in the range from 10 kHz to about 200 kHz. In certain implementations, when the systems and methods are configured for radar imaging, the frequencies may be in the range from 100 MHz to about 30 GHz. Generally, the systems and methods described herein may be applied for any frequency range, without departing from the scope of the present disclosure.

Overpinging with Multiple SAS Transmitters

In water, the maximum range of a sonar is typically defined as the distance sound can travel to a target and back before the next transmission. However, after the next transmission, the prior ping continues to propagate through the water. If it is possible to use a prior ping, then the range of a sonar and/or area coverage rate of a sonar may be increased. If N total pings (N−1 prior pings) are fully used, then the area coverage rate may be increased by a factor of N.

Unfortunately, distant echoes are weaker than closer ones (assuming constant target strength). By using orthogonal signals, it may be possible to improve the signal to noise ratio (SNR), but it will still be nearly impossible to receive while a transmitter is transmitting. For a system with a classical range of R, this means it may be (almost) always impossible to observe echoes at ranges just past R.

Accordingly, in yet another aspect, the present disclosure relates to a solution to this SNR problem by adding multiple transmitters that are spaced and fired in a novel order so as provide multiple opportunities to recover data from the same range and phase center position and ensure that at least one of those observations is not jammed. In some implemenations, N separated transmitters may be used to increase the range of a virtual array by N.

In one example, consider an array with length L and maximum classical SAS range R at a velocity V (V*dt=L/2, R=c/2*dt=cL/4V) with a delay between pings of dt=L/2V. Assume that transmitters are (with zero being on the bow and distance increasing moving aft) at x=0, x=0.25L, x=0.6L. For the second transmitter to form the same virtual array as the first it must delay firing until the vehicle has moved sufficiently to position it (corresponding to about 0.25 dt). Likewise, if the third transmitter is used, it must be delayed 0.6 dt. Suppose ping 1 is formed using transmitters 1 and 2, ping 2 is formed using transmitters 1 and 3, ping three is formed using transmitters 2 and 3, and then the sequence is repeated. The firing sequence timing may then be:

Transmitter 1: [0, 1, off, 3, 4, off, . . . ]*dt

Transmitter 2: [0.25, off, 2.25, 3.25, off, 5.25, . . . ]*dt

Transmitter 3: [off, 1.6, 2.6, off, 4.6, 5.6, . . . ]*dt

The combination of all ping times is therefore: [0, 0.25, 1, 1.6, 2.25, 2.6, 3, 3.25, 4, 4.6, 5.25, 5.6, etc] *dt. Transmitter 1 and 2 each form virtual arrays for ping position 1. Since jamming is caused by future pings, virtual array 1 ping 1 is jammed at the following ranges: [0.25, 1, 1.6, 2.25, 2.6, 3, 3.25, 4, 4.6, 5.25, 5.6, etc]*R. Virtual array 2 ping 1 is jammed at the following ranges: ([1, 1.6, 2.25, 2.6, 3, 3.25, 4, 4.6, 5.25, 5.6, etc]−0.25)*dt=[0.75, 1.35, 2, 2.35, 2.75, 3, 3.75, 4.35, 5, 5.35, etc] *R.

Since virtual array 1 is jammed at 0.25 R, but virtual array 2 is not, virtual array 2's signal is used for those immediate ranges. Since virtual array 2 is jammed at 0.75 R but virtual array 1 is not, virtual array 1's signal is used there. Since both arrays are jammed at 2R it is not possible to get an unjammed observation of 2R, therefore that is the greatest unjammed range. The ranges for which both arrays are jammed (the ranges at which the fused signal is jammed) are: [3,6,9,12 etc] *R. Likewise, the second ping uses virtual arrays formed by transmitters 1 and 3, the third ping uses arrays 2 and 3, etc. Signals are interlaced in a similar matter for latter pings as well. Changing the signals from ping to ping further reduces noise suppression.

Although the transducers/projectors spacing was described above as being 0, 0.25*L and 0.6*L, the system 100 may include any number of projectors positioned at any suitable spacing and having any suitable firing sequences without departing from the scope of the present disclosure. Generally, system 100 may combine overpinging with multiple alongtrack projectors by either repeating the patter or by generating a new sequence, which may be random. As another example, system 100 may include an array of five projectors positioned at 0, 0.25*L, 0.4142*L, 0.6*L and 0.7321*L. Position 0.4142*L corresponds to [sqrt(2)−1]*L and position 0.732*L corresponds to [sqrt(3)−1]*L. In such an example, the firing sequences may be:

Transmitter 1: [0, 1, 2, 3, off, 5 etc.] *dt

Transmitter 2: [0.25, 1.25, 2.25, off, 4.25, 5.25, etc.] *dt

Transmitter 3: [0.4142, 1.4142, off, 3.4142, 4.4142, 5.4142, etc.] *dt

Transmitter 4: [0.6, off, 2.6, 3.6, 4.6, 5.6, etc.] *dt

Transmitter 5: [off, 1.7321, 2.7321, 3.7321, 4.7321, off, etc.] *dt

They systems and methods described herein may be adapted as desired for both sonar and radar systems, and accordingly for both synthetic aperture sonar (SAS) and synthetic aperture radar (SAR) systems. For example, sonar transducers may be replaced with suitable radar transducers, and one or more components may be modified, added to or removed from the systems described herein to operate in a sonar and radar regime. In some implementations, the systems and methods may be configured to operate as both sonar and radar devices, without departing from the scope of the present disclosure. In certain implementations, when the systems and methods are configured for sonar imaging, the frequencies may be in both high and low frequency ranges in the range from 10 kHz to about 200 kHz. In certain implementations, when the systems and methods are configured for radar imaging, the frequencies may be in the range from 100 MHz to about 30 GHz. Generally, the systems and methods described herein may be applied for any frequency range, without departing from the scope of the present disclosure.

Holographic Simultaneous Localization and Mapping (SLAM)

Classical holographic navigation for AUVs, by requiring at least one image to be from a synthetic aperture, does not enable a true SLAM (simultaneous localization and mapping) solution. This is because when a real aperture image is correlated with a synthetic aperture image, the position estimate update is an average one; the estimate can not be isolated to individual states.

In one aspect, the present disclosure relates to a method of removing the above described requirement of needing at least one synthetic aperture image or quasi-hologram for holographic navigation. In some implementations, region connecting observation positions are defined and used in such a way as to enable real aperture correlation, e.g., by forming “correlation tubes.” An idealized sonar is introduced and improves the correlation performance inside the tubes. In some implementations, a new method of operation allows a typical survey sonar to use correlation tubes to improve its navigation.

As noted earlier, generally, holographic navigation works because a hologram which is defined over a range of angles contains all possible images within that set of angles (subject to a few constraints such as frequency limitations, occlusion, etc). When a real aperture image is correlated against a synthetic aperture image, the correlation process transparently identifies the position in the quasi-hologram where the real aperture image originated. The correlation works because although the real aperture sonar sees each object from one vantage point, and sees a set of objects from a set of different vantage points, the quasi-hologram is a multi-aspect record. When two real aperture images are compared, only targets along a line connecting the two sonars can be correlated, previously assumed to be impossible. Since the percentage of features that are collinear with the two sonars is exceptionally small, and since the correlation result from the area off of that line is noise, the resulting signal to noise ratio for the correlation is near zero.

In some implementations, the above described method solves this problem by combining a preliminary navigation estimate, a real aperture array, and “correlation tubes.” In some implementations, if the AUV has a reasonably accurate estimate of its position when it made a prior observation, it can define a reasonably accurate vector connecting the two positions. Using a real aperture array, the sonar can steer the signals received at either position along the direction of that vector. By forming both beams down the correct correlation direction, the signal to noise ratio improves considerably.

In some implementations, the real aperture sonars are long enough that both are in the nearfield and it is possible to form a beam without spurious information (there is a temptation to call it “noise free”, but it would, of course, still be subject to environmental noise sources). For example, consider two parallel arrays of length L separated by a distance D along the broadside vector of the arrays. If the resolution of the arrays at the separation distance is less than the length of the arrays then the arrays have sufficient resolution to “block out” energy from areas with non-overlapping aspects. (Mathematically, this is very similar to a nearfield constraint). Assuming a wavelength λ and making a small angle assumption, the angular resolution of the arrays is ΔΘ=2/L, and the across range resolution is Δx=rΔΘ=r λ/L. Since we want Δx<L or r λ/L<L, the technique performs best when r<L²/λ. The nonlinear relationship between r and L means that the maximum range increases quickly as the array length grows. For most existing synthetic aperture sonars, with L≈0.50 m and λ≈0.01 m, the maximum range is approximately 25 m. In some implementations, with L=2.5 m and λ=0.0083 m the maximum range is approximately 753 m, or three times the intended survey range of 250 m. More generally, given a beam width ϕ the range constraint becomes r<cos² ϕL²/λ. According to an illustrative implementation, a system designed according to the beam width constraint can use correlation tubes off of broadside without performance degradation.

FIG. 12 depicts a process for simultaneous localization and mapping (SLAM) using real aperture sonar images, according to an illustrative implementation of the present disclosure. The sonar system 100 on a vehicle may receive a first real aperture acoustic data of a portion of the terrain being traversed by the vehicle (step 1202). The sonar system 100 (e.g., CCU 160) may receive a first position estimate representing a position from which the first real aperture image was obtained (step 1204). In some instances, the position estimate and the first data (or image) may have been obtained a priori (either on the same mission or on a prior mission) by the same vehicle as it traverses the terrain. In other instances, the position estimate and the first data (or image) may have been obtained by another vehicle (at the same time or different time) traversing the terrain. Both vehicles may be traversing the terrain simultaneously and communicating with each other. The vehicle receiving the images may be moved to a current position (step 1206). At the current position, the vehicle may receive a second position estimate representing the current position (step 1208). The sonar system 100 may then determine the correlation axis (or correlation tube) connecting the first position estimate and the second position estimate (step 1210). The correlation may determine the direction of beamforming for the vehicle in the current position. The sonar system 100 may insonify the terrain and generate a second real aperture image of the terrain from the current position (step 1212). The beamformer 134 of the sonar system 100 may steer the receiving signal such that the received acoustic signal is directed towards the terrain and along the correlation axis or tube (step 1212), thereby allowing the vehicle at the current position to view the terrain along the same axis as the direction along which the first real aperture image was obtained.

Most synthetic aperture sonars are used to perform surveys. Most surveys follow a lawnmower pattern. Although it is typical when surveying using a lawnmower pattern to overlap adjacent passes (to ensure full coverage), that overlap is not suitable for holographic navigation (since it has the wrong spectral orientation). Assuming the sonar has a maximum useful range R, swaths are typically spaced 2R apart, and the portion of the imagery from an adjacent survey leg that is oriented correctly for holographic navigation exists at the ranges between 2R and 3R. The sonar design of the systems herein is advantageous as it enables clean correlation at triple the baseline operating range.

In some implementations, there may be a plurality of ways to navigate relative to distant sonar imagery (without looping around to look at it). One such way may include a SAS, such as sonar system 100 and slowing the vehicle down to ⅓ speed, which will then triple the range, enabling a holographic navigation fix. Generally, the vehicle may be slowed down any amount suitable to increase the range by a desirable amount. Another way is to include a dual frequency system, one frequency can map continuously while another alternates between a normal ping rate (to map) and a ⅓ ping rate (to observe distant regions). A third way may be to use sequences of orthogonal signals. If the distant echoes are sufficiently orthogonal to the closer echoes they can be used for correlation and a navigation update. Unlike imaging at long range, it is not necessarily to use signals taken at all ranges (including portions where the SNR is poor due to other transmissions). It is generally sufficient to use only range regions with “good enough” SNR. When correlating relative to a prior pass it is possible to steer beams through many different prior robot positions, producing many different measurements. Used in post-processing this can result in an accurate map, although, in some implementations, a limited set of beams may be processed in real time.

They systems and methods described herein may be adapted as desired for both sonar and radar systems, and accordingly for both synthetic aperture sonar (SAS) and synthetic aperture radar (SAR) systems. For example, sonar transducers may be replaced with suitable radar transducers, and one or more components may be modified, added to or removed from the systems described herein to operate in a sonar and radar regime. In some implementations, the systems and methods may be configured to operate as both sonar and radar devices, without departing from the scope of the present disclosure. In certain implementations, when the systems and methods are configured for sonar imaging, the frequencies may be in both high and low frequency ranges in the range from 10 kHz to about 200 kHz. In certain implementations, when the systems and methods are configured for radar imaging, the frequencies may be in the range from 100 MHz to about 30 GHz. Generally, the systems and methods described herein may be applied for any frequency range, without departing from the scope of the present disclosure.

Bistatic and Monostatic Gapfiling for SAS

Synthetic aperture sonars perform poorly in the near “nadir” regime (directly under the vehicle). This is because the partial derivative of range with respect to horizontal distance may be approximately zero directly beneath the sonar.

There are two traditional solutions to the aforementioned problem. The first is to use a real aperture sonar to image directly under the vehicle. However, as the survey sonar range increases it is generally necessary to survey from a higher altitude, causing the resolution of a real aperture gapfiller generally to drop for a given aperture. This is due to two effects: decreased spatial resolution due to increased range, and needing to use lower frequencies/longer wavelengths due to range and absorption. Even though the aperture can be increased, at realistic altitudes for a long range survey sonar, it is impossible to get resolutions that approach SAS resolutions.

Applicants solve this problem by using a plurality of sonar vehicles that are positioned so as to be able to hear each other's transmissions and image bistatically. The bistatic image has SAS level resolution under either vehicle, and poor resolution in the gaps between the vehicles. By fusing SAS imagery and bistatic SAS imagery it is possible to have a very high resolution map.

In certain implementations, a long range vehicle may be used in conjunction with a smaller vehicle. In such implementations, the long range vehicle used for imaging, may include a very large gap beneath it. A second, smaller vehicle may then be used specifically to image the large vehicle's gap. The second vehicle may fly the same track line, but below the large vehicle.

In certain other implementations, a vehicle may be used at higher altitude and create gap beneath the vehicle that is comparable to the range of the sonar. In such implementations, one or more adjacent mission legs are then flown such that subsequent legs fill the gap of prior legs with minimal waste. This design assumes a steeper grazing angle than traditional SAS, but would be considered more typical of synthetic aperture radar (SAR). From a grazing angle compensation/holographic navigation perspective, described above, the reduced range of grazing angles may be beneficial.

Passive Relative Localization

When multiple vehicles perform a survey it is necessary to overlap their survey areas so as to account for inertial navigation errors that accrue over the mission. (This assumes an unknown area without a predeployed beacon network such as long baseline or ultrashort baseline.) Since the inertial navigation error grows with time, there can be substantial drift over a long mission. This large drift requires a large overlap, substantially decreasing the net area coverage rate of the sonar. If, instead, vehicles can fly in tight formation then that overlap can be reduced, and gross errors only occur at the edge of the areas imaged by the formation.

One way for the vehicles to maintain their formation is by using beacon systems. Using an onboard ultrashort baseline navigation system, it is possible to measure the range and bearing to another vehicle. However, this requires an additional system.

Instead, Inventor proposes a reduction in which the vehicle (which is assumed to have a synthetic aperture sonar array and system 100) passively listens for the transmissions from other vehicles using its real aperture sonar. Using the received signal it may measure the bearing to the other vehicle (but not the range). Measuring range passively is difficult since most SAS's are “slaved to speed” meaning that the pings are timed based on the perceived position of the vehicle. Schemes may be used to passively estimate range based on time of arrival, but the random component of SAS ping timing makes this undesirable.

The methods and systems described herein measure range by passively dithering the vehicles relative to one another and fusing the data in a navigation filter. For instance, consider two vehicles flying in parallel with slowly drifting inertial navigation systems. Assume that they each have a base survey velocity of 2 m/s. Assume that for 5 minutes vehicle #1 flies at 2.05 m/s and vehicle #2 flies at 1.95 m/s, resulting in a 30 m change in position. If the vehicles are 300 m apart, this corresponds to a 5.6 degree change in position; if they are 310 m apart, this may correspond to an 5.8 degree change. A 2.5 m array with a wavelength of 0.008 m has an angular resolution of 0.18 degrees, making a 10 m range variation observable.

In an alternative instantiation of methods and systems described herein, the vehicle passively listens to a timed pinger on other vehicles to estimate range. A filter onboard the vehicle estimates clock drift. Although dithering is not necessary if listening to a passive pinger dithering does make clock drift more observable.

Bi-Static and Multi-Static Holographic Navigation

In another aspect, the systems and methods described herein relate to localizing an emitter or receiver with high precision relative to the sea floor. It can be used as either a fully active system (using well synchronized transmitters and receivers) or as a passive system (localizing a transmitter or a receiver based on poorly timed received signals).

The system has many applications. For instance, consider the case of localizing a vehicle based on transmissions from a known source. This applies to scenarios such as mine neutralization and submarine navigation. Bistatic holographic navigation allows small autonomous mine neutralizers to use low frequency signals. Low frequency acoustic signals are more robust to sedimentation than high frequency signals, but low frequency projectors can weigh more than neutralization platforms.

Placing the low frequency projector onto another platform (such as a large diameter unmanned underwater vehicle) allows the mine neutralizer to reduce its sonar to small number of receive elements (which weigh very little). This enables low frequency holographic navigation on a very small platform and, ultimately, a more robust and flexible mine neutralization solution. Similarly, submarines prefer not to use active sonar as it reveals their location, but they have very high quality receivers. Bi-static holographic navigation allows submarines to localize themselves based on seafloor scattering from a known source (either fixed or moving).

Bi-static holographic navigation can also be used to localize a transmitter using a known receiver. Submarines and stealth aircraft often possess sonars and radars with electronically scanned phased arrays that steer energy into narrowly defined areas to reduce the probability of intercept. Although their outbound beams are narrow, the scattering off the seafloor or ground is less directional. Bi-static holographic navigation can detect and localize those target based on scattering. Similarly, ships and whales are difficult to accurately track and localize; bistatic holographic navigation can be used to track them and reduce collisions (a half dozen blue whales were killed off of San Francisco last year).

Bi-static holographic navigation, in its simplest form, may include two platforms with accurate synchronized clocks, a known transmit signal with known timing, position information about the transmitter, and a prior bi-static map. The receive system (either a single channel or an array) receives the signal, forms an image using timing and geometry information, applies grazing angle compensation, correlates the image against a prior image, converts the correlation result to a probability density function, and uses that probability density function to update its estimate of the state (position). HF holographic navigation techniques can be applied to remove various forms of phase error. In a more advanced form, the emitting system (submarine, aircraft, etc) would emit a known signal from an unknown position at an unknown time that would be received at a known location at a known time. That signal would be back-projected or beamformed, correlated, and used to update a position estimate.

In a more difficult case, the signal being emitted by the submarine or aircraft would be unknown. In this case, if the direct arrival were known, it could be used to form an image. If the signal were unknown, it would be necessary to learn the signal, use it to form an image, and localize the target.

These techniques can also be performed multi-statically and, as always, apply equally to sonar and radar. Although back-scatter is more desirable, forward scatter can be used. Bi-static holographic navigation follows a similar process to monostatic holographic navigation. First, it is necessary to create a prior map. Ideally, this is a bi-static synthetic aperture image containing all possible frequencies at all possible angles. In practice, this will likely be band limited and angle limited. The remaining process depends on whether the transmitter and receiver are synchronized, which system has a known position, and whether the signal is known. It is possible, although undesirable, to use two unknown positions.

FIG. 13 shows a bi-static holographic system 1300 including a transmitter 1302 and receivers 1304 and 1314 operating in an underwater environment 1318. The transmitter 1302 may reside on a first platform such as an underwater vessel or surface vessel. The receivers 1304 and 1314 may reside on second and third platforms that include an AUV, submarine, robot, and the like. FIG. 13 illustrates an acoustic signal 1310 directed from transmitter 1302 that results in a scattering signal 1312 reflected from a portion 1316 of the underwater terrain 1308. A portion 1312 a of the scattering signal is directed toward the receiver 1304 while another portion 1312 c is directed toward receiver 1314. While FIG. 13 illustrates a system 1300 for navigating an underwater terrain 1308, other systems may utilize the similar technique for navigation using radar instead of sonar.

The acoustic transmitter 1302 may be stationary or in motion. The receiver 1304 or 1314 may be stationary or in motion. The position of the transmitter 1302 may be known by the receivers 1304 and/or 1314. The position of the receivers 1304 and/or 1314 may be known by the transmitter 1302. While FIG. 1 refers to a system 100 capable of performing both transmit and receiver functions associated with holographic navigation, one of ordinary skill will readily recognize certain functions of the system 100 can be implemented on separate a transmitter 1302 and receiver 1304 platforms to enable bi-static or multi-static holographic navigation. For example, the transmitter 1302 may include portions of the sonar unit 110 such as transducer array 112, transmitter 116, and transmitter controller 118. The transmitter 1302 may include a central controller 160 or processor that controls operations of the transmitter 1302. The transmitter may include a memory that store certain data such as data in map store 154. The transmitter 1302 may include a clock, or the central control unit 160 or navigation controller 170 may implement a clock. The transmitter 1302 may include any other components of system 100 to enable operation of the transmitter 1302, or of the platform on which the transmitter resides.

One of ordinary skill will also recognize that certain functions of the system 100 can be implemented on a separate receiver 1304 or 1314 to enable bi-static or multi-static navigation. For example, the receiver 1304 or 1314 may include portions of the sonar unit 110 such as transducer array 112 and receiver 114. The receiver 1304 or 1314 may include a preprocessor 120 including filter conditioner 122 and Doppler compensator 124. The receiver 1304 or 1314 may include a matched filter 130 having a pulse compensator 132 and beamformer 134. The receiver 1304 or 1314 may include a signal corrector 104 having a grazing angle compensator 142 and phase error corrector 144. The receiver 1304 or 1314 may include a signal detector 150 having a signal correlator 152 and map store 154 or other memory to store navigation-related data. The receiver 1304 or 1314 may include a clock, or the central control unit 160 or navigation controller 170 may implement a clock.

The receiver 1304 or 1314 may include any other components of system 100 to enable operation of the receiver 1304 or 1314, or of the platform on which the receiver 1304 or 1314 resides. The clock of the transmitter 1302 may be synchronized with a clock at the receiver 1304 and/or 1314. Synchronization may be performed prior to launch of the transmitter and receiver platforms or during operation of the platforms. The platform may communication information with each other to enable clock synchronization, update location information, and/or to exchange map data. Communications between platforms may be implemented via direct acoustic signals, direct wireless link while the platforms are surfaced, or indirect wireless links (e.g., satellite links, ship-to-ship links, and the like).

Depending on whether the transmitter 1302 and receiver 1304 or 1314 are synchronized, which platform (transmitter or receiver) has a known position, and whether the signal is known, the system 1300 is able to estimate the position of either the transmitter 1302 or receivers 1304 and 1314. The system 1300 can also estimate the position of either the transmitter 1302 or receivers 1304 and 1314 when the positions of the transmitter 1302 or receivers 1304 and 1314 is unknown. FIGS. 14-19 are flow diagrams that illustrate the how the system 1300 estimates the positions of the receiver 1304 or 1314, or transmitter 1302 depending on the type of information known by the receiver 1304 or 1314 and/or the transmitter 1302. An acoustic signal may include one or more characteristics or settings such as, without limitation, a noise sequence, linear chirp, encoding, power, frequency, phase, and the like.

FIG. 14 is a flow diagram of a process 1400 for estimating a receiver 1304 or 1314 position where timing is known, transmitter position is known, the acoustic signal 1310 is known, but the receiver 1304 or 1314 position is unknown. In this scenario, the receiver 1304 or 1314 has its clock synchronized with a clock of the transmitter 1302 so that the timing of the acoustic signal is known by the receiver 1304 or 1314. The transmitter 1302 transmits an acoustic signal 1310 through the underwater environment 1318 toward the portion 1316 of the underwater terrain 1308. The acoustic signal 1310 is reflected or scattered from the portion 1316 of the underwater terrain 1308 in multiple directions. A portion 1312 a of the scattered signal is received by the receiver 1304 (block 1402). Again, the receiver may be an AUV or submarine. The receiver 1304 applies the received signal 1312 a to a matched filter 130 (block 1404). The receiver 1304 may also apply the signal to a preprocessor 120 prior to the matched filter 130. The receiver forms an image using, for example, beamformer 134 (block 1406).

Once an image is formed, the receiver applies bi-static grazing angle compensation using grazing angle compensator 142 (block 1408). The image may also be phase error corrected by a phase error corrector 144 when high frequency holographic navigation is used. The receiver 1304 then correlates the image with a prior map as discussed previously herein with respect to holographic navigation (1410). The map may include a monostatic or bi-static map. Over modest angles, a monostatic prior map can be a substitute for a bistatic prior map. The receiver 1304 may use a processor such as central control unit 160 or a dedicated processor to convert the correlation result to a probability density function (block 1412). The receiver 1304, using a processor such as central control unit 160, then updates the state estimate (e.g., location) of receiver 1304 (block 1414). The processor 160 may also determine the heading or direction of motion of the receiver 1304 and/or platform on which the receiver 1304 resides.

In certain implementations, the real image is taken and a sum of the envelopes of small image correlation regions with approximately stationary phase is calculated before calculating a probability density function based on that sum. This may be similar to a speckle reduction technique used in imaging methods. The sums can either be for a single altitude solution or for multiple altitude solutions; if multiple altitude solutions are used then the technique measures altitude bias. In certain implementations, using the envelope only (the absolute value of the correlation result) removes the relative phase differences between correlation results. It is important to note that summing together a large number of correlation images results in a transition from Rayleigh to Gaussian distributed speckle intensity. The receiver 1304 may account for this difference when converting the correlation result to a probability density function. Previous holographic navigation techniques, which use the Rayleigh distribution, typically fail when presented with correlation results based on sub-image summation. By switching to a more representative distribution, the receiver 1304 avoids this deficiency. When a small number of images are summed together the distribution may not yet be fully Gaussian. Hence, the receiver 1304 may be better represented the distribution using some other distribution such as a K-distribution.

FIG. 15 is a flow diagram of a process 1500 for estimating a receiver 1304 or 1314 position where timing is unknown, transmitter 1302 position is known, the acoustic signal 1310 is known and the receiver 1304 or 1314 position is unknown. In this scenario, the timing of the acoustic signal is unknown by the receiver 1304 or 1314. The transmitter 1302 transmits an acoustic signal 1310 through the underwater environment 1318 toward the portion 1316 of the underwater terrain 1308. The acoustic signal 1310 is reflected or scattered from the portion 1316 of the underwater terrain 1308 in multiple directions. A portion 1312 a of the scattered signal is received by the receiver 1304 (block 1502). The receiver 1304 applies the received signal 1312 a to a matched filter 130 (block 1504). The receiver 1304 may also apply the signal to a preprocessor 120 prior to the matched filter 130. The receiver forms multiple images for hypothesis testing using, for example, beamformer 134 (block 1506).

Once the multiple images are formed, the receiver applies bi-static grazing angle compensation using grazing angle compensator 142 (block 1508). The images may also be phase error corrected by a phase error corrector 144 when high frequency holographic navigation is used. The receiver 1304 then correlates the images with a prior map as discussed previously herein with respect to holographic navigation (block 1510). The receiver 1304 may use a processor such as central control unit 160 or a dedicated processor to convert the correlation result to a probability density function (block 1512). The receiver 1304, using a processor such as central control unit 160, determines which of the multiple images is most accurate based on hypothesis testing and then updates the state estimate (e.g., location) of receiver 1304 (block 1514). Hypothesis testing may include a processing in which the receiver 1304 makes crude estimate of the transmitter 1302 location and then creates different images based on where the transmitter could be located. The processor 160 may also determine the heading or direction of motion of the receiver 1304 and/or platform on which the receiver 1304 resides.

FIG. 16 is a flow diagram of a process 1600 for estimating a receiver 1304 position where timing is unknown, transmitter 1302 position is known, the acoustic signal 1310 is unknown and the receiver 1304 position is unknown. In this scenario, the transmitter 1302 transmits an acoustic signal 1310 through the underwater environment 1318 toward the portion 1316 of the underwater terrain 1308. The acoustic signal 1310 is reflected or scattered from the portion 1316 of the underwater terrain 1308 in multiple directions. A portion 1312 a of the scattered signal is received by the receiver 1304 (block 1602). The receiver 1304 learns the acoustic signal 1310 based on detecting one or more of the portions of the scattered signal 1312 and other known data such as the transmitter 1302 position (block 1604). The receiver 1304 may employ any number of learning techniques such as, without limitation, a technique including a simulated annealing algorithm and/or process to learn the acoustic signal 1310. The receiver 1304 applies the received signal 1312 a to a matched filter 130 (block 1606). The receiver 1304 may also apply the signal to a preprocessor 120 prior to the matched filter 130. The receiver forms an image using, for example, beamformer 134 (block 1608).

Once an image is formed, the receiver applies bi-static grazing angle compensation using grazing angle compensator 142 (block 1610). The image may also be phase error corrected by a phase error corrector 144 when high frequency holographic navigation is used. The receiver 1304 then correlates the image with a prior map as discussed previously herein with respect to holographic navigation (1612). The receiver 1304 may use a processor such as central control unit 160 or a dedicated processor to convert the correlation result to a probability density function (block 1614). The receiver 1304, using a processor such as central control unit 160, then updates the state estimate (e.g., location) of receiver 1304 (block 1616). The processor 160 may also determine the heading or direction of motion of the receiver 1304 and/or platform on which the receiver 1304 resides.

FIG. 17 is a flow diagram of a process 1700 for estimating a transmitter 1302 position where timing is known, transmitter 1302 position is unknown, the acoustic signal 1310 is known and the receiver 1304 position is known. In this scenario, the receiver 1304 or 1314 has its clock synchronized with a clock of the transmitter 1302 so that the timing of the acoustic signal 1310 is known by the receiver 1304 or 1314. The transmitter 1302 transmits an acoustic signal 1310 through the underwater environment 1318 toward the portion 1316 of the underwater terrain 1308. The acoustic signal 1310 is reflected or scattered from the portion 1316 of the underwater terrain 1308 in multiple directions. A portion 1312 a of the scattered signal is received by the receiver 1304 (block 1702). Again, the receiver may be an AUV or submarine. The receiver 1304 applies the received signal 1312 a to a matched filter 130 (block 1704). The receiver 1304 may also apply the signal to a preprocessor 120 prior to the matched filter 130. The receiver forms an image using, for example, beamformer 134 (block 1706).

Once an image is formed, the receiver applies bi-static grazing angle compensation using grazing angle compensator 142 (block 1708). The image may also be phase error corrected by a phase error corrector 144 when high frequency holographic navigation is used. The receiver 1304 then correlates the image with a prior map as discussed previously herein with respect to holographic navigation (block 1710). The receiver 1304 may use a processor such as central control unit 160 or a dedicated processor to convert the correlation result to a probability density function (block 1412). The receiver 1304, using a processor such as central control unit 160, then updates the state estimate (e.g., location) of transmitter 1302 (block 1714). The processor 160 may also determine the heading or direction of motion of the transmitter 1302 and/or platform on which the transmitter 1302 resides.

FIG. 18 is a flow diagram of a process 1800 for estimating a transmitter 1302 position where timing is unknown, transmitter 1302 position is unknown, the acoustic signal 1301 is known and the receiver 1314 position is known. In this scenario, the timing of the acoustic signal 1310 is unknown by the receiver 1304 or 1314. The transmitter 1302 transmits an acoustic signal 1310 through the underwater environment 1318 toward the portion 1316 of the underwater terrain 1308. The acoustic signal 1310 is reflected or scattered from the portion 1316 of the underwater terrain 1308 in multiple directions. A portion 1312 a of the scattered signal is received by the receiver 1304 (block 1802). The receiver 1304 applies the received signal 1312 a to a matched filter 130 (block 1804). The receiver 1304 may also apply the signal to a preprocessor 120 prior to the matched filter 130. The receiver forms multiple images for hypothesis testing using, for example, beamformer 134 (block 1806).

Once the multiple images are formed, the receiver applies bi-static grazing angle compensation using grazing angle compensator 142 (block 1808). The images may also be phase error corrected by a phase error corrector 144 when high frequency holographic navigation is used. The receiver 1304 then correlates the images with a prior map as discussed previously herein with respect to holographic navigation (block 1810). The receiver 1304 may use a processor such as central control unit 160 or a dedicated processor to convert the correlation result to a probability density function (block 1812). The receiver 1304, using a processor such as central control unit 160, determines which of the multiple images is most accurate based on hypothesis testing and then updates the state estimate (e.g., location) of transmitter 1302 (block 1814). The processor 160 may also determine the heading or direction of motion of the transmitter 1302 and/or platform on which the transmitter 1302 resides.

FIG. 19 is a flow diagram of a process 1900 for estimating a transmitter 1302 position where timing is unknown, transmitter 1302 position is unknown, the acoustic signal 1310 is unknown, but the receiver position 1304 or 1314 is known. In this scenario, the transmitter 1302 transmits an acoustic signal 1310 through the underwater environment 1318 toward the portion 1316 of the underwater terrain 1308. The acoustic signal 1310 is reflected or scattered from the portion 1316 of the underwater terrain 1308 in multiple directions. A portion 1312 a of the scattered signal is received by the receiver 1304 (block 1902). The receiver 1304 learns the acoustic signal 1310 based on detecting one or more of the portions of the scattered signal 1312 and other known data such as the transmitter 1302 position (block 1904). The receiver 1304 applies the received signal 1312 a to a matched filter 130 (block 1906). The receiver 1304 may also apply the signal to a preprocessor 120 prior to the matched filter 130. The receiver forms an image using, for example, beamformer 134 (block 1908).

Once an image is formed, the receiver 1304 applies bi-static grazing angle compensation using grazing angle compensator 142 (block 1910). The image may also be phase error corrected by a phase error corrector 144 when high frequency holographic navigation is used. The receiver 1304 then correlates the image with a prior map as discussed previously herein with respect to holographic navigation (block 1912). The receiver 1304 may use a processor such as central control unit 160 or a dedicated processor to convert the correlation result to a probability density function (block 1914). The receiver 1304, using a processor such as central control unit 160, then updates the state estimate (e.g., location) of transmitter 1302 (block 1916). The processor 160 may also determine the heading or direction of motion of the transmitter 1302 and/or platform on which the transmitter 1302 resides.

The systems and methods described herein may be realized as a software component operating on a conventional data processing system such as a Unix system. In that implementation, these mechanisms can be implemented as a C language computer program, or a computer program written in any high level language including Matlab, C++, Fortran, Java or BASIC. Additionally, in an implementation where microcontrollers or DSPs are employed, the mapping mechanism can be realized as a computer program written in microcode or written in a high level language and compiled down to microcode that can be executed on the platform employed. The development of such data processing systems is known to those of skill in the art, and such techniques are set forth in Digital Signal Processing Applications with the TMS320 Family, Volumes I, II, and III, Texas Instruments (1990). Additionally, general techniques for high level programming are known, and set forth in, for example, Stephen G. Kochan, Programming in C, Hayden Publishing (1983). It is noted that DSPs are particularly suited for implementing signal processing functions, including preprocessing functions such as image enhancement through adjustments in contrast, edge definition and brightness. Developing code for the DSP and microcontroller systems follows from principles well known in the art. The system also provides and enables as is known to those of skill in the art, object oriented frameworks are generally understood as a set of classes that embody a design for solutions to a family of related problems. See The C++ Programming Language, 2nd Ed., Stroustrup Addision-Wesley. Accordingly, a framework for mapping and filtering may be created that provides a prefabricated structure, or template, of a working mapping and filtering program.

Variations, modifications, and other implementations of what is described may be employed without departing from the spirit and scope of the systems and methods described herein. For example, though the systems and methods are described in the context of underwater mapping and navigation using sonar signals, the systems and methods may be equally applicable for mapping and navigating in aerial or other land or space-based terrains and using other imaging technologies include radar, optical signals, and any acoustic or electromagnetic signal. Moreover, any of the method and system features described above or incorporated by reference may be combined with any other suitable method or system feature disclosed herein or incorporated by reference, and is within the scope of the contemplated systems and methods. The systems and methods may be embodied in other specific forms without departing from the spirit or essential characteristics thereof. The foregoing implementations are therefore to be considered in all respects illustrative, rather than limiting of the systems and methods described herein. 

What is claimed is:
 1. A high frequency sonar system comprising a phase error corrector configured to: receive a sonar image; estimate a range varying phase error of the received sonar image; determine that a variation in the range varying phase error across the received sonar image is greater than a first error threshold; in response to determining that the variation in the range varying phase error across the received sonar image is greater than the first error threshold: split the received sonar image into a plurality of sub-regions; estimate, for each sub-region of the plurality of sub-regions, a respective range varying phase error; determine, for a sub-region of the plurality of sub-regions, that a variation in the range varying phase error across a respective sub-region is less than a second error threshold; in response to determining that the variation in the range varying phase error across the respective sub-region is less than the second error threshold, correct the sub-region to reduce the variation in the range varying phase error across the respective sub-region.
 2. The system of claim 1, wherein the received sonar image is a real aperture image.
 3. The system of claim 2, wherein the phase error corrector is further configured for modifying the real aperture image with grazing angle compensation.
 4. The system of claim 2, wherein the first error threshold is based on a maximum range of the real aperture image.
 5. The system of claim 1, wherein the received sonar image is a synthetic aperture image.
 6. The system of claim 1, wherein the sub-regions are of equal size.
 7. The system of claim 1, wherein the sub-regions are of varying size, and the variation in size is based on the range.
 8. The system of claim 1, wherein the first error threshold is the same as the second error threshold.
 9. The system of claim 1, wherein the sub-regions have constant altitude phase error or constant sound speed error.
 10. The system of claim 1, wherein the respective range varying phase error is substantially constant across the sub-region.
 11. A method of correcting range varying phase errors, comprising: receiving a sonar image; estimating a range varying phase error of the received sonar image; determining that a variation in the range varying phase error across the received sonar image is greater than a first error threshold; in response to determining that the variation in the range varying phase error across the received sonar image is greater than the first error threshold: splitting the received sonar image into a plurality of sub-regions; estimating, for each sub-region of the plurality of sub-regions, a respective range varying phase error; determining, for a sub-region of the plurality of sub-regions, that a variation in the range varying phase error across a respective sub-region is less than a second error threshold; in response to determining that the variation in the range varying phase error across the respective sub-region is less than the second error threshold, correcting the sub-region to reduce the variation in the range varying phase error across the respective sub-region.
 12. The method of claim 11, wherein the received sonar image is a real aperture image.
 13. The method of claim 12, further comprising modifying the real aperture image with grazing angle compensation.
 14. The method of claim 12, wherein the first error threshold is based on a maximum range of the real aperture image.
 15. The method of claim 11, wherein the received sonar image is a synthetic aperture image.
 16. The method of claim 11, wherein the sub-regions are of equal size.
 17. The method of claim 11, wherein the sub-regions are of varying size, and the variation in size is based on the range.
 18. The method of claim 11, wherein the first error threshold is the same as the second error threshold.
 19. The method of claim 11, wherein the sub-regions have constant altitude phase error or constant sound speed error.
 20. The method of claim 11, wherein the respective range varying phase error is substantially constant across the sub-region. 